Nonlinear gravitational self-force. I. Field outside a small body 
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A small extended body moving through an external spacetime g a p creates a metric perturbation 
h a p, which forces the body away from geodesic motion in g a p. The foundations of this effect, 
called the gravitational self-force, are now well established, but concrete results have mostly been 
limited to linear order. Accurately modeling the dynamics of compact binaries requires proceeding 
to nonlinear orders. To that end, I show how to obtain the metric perturbation outside the body 
at all orders in a class of generalized wave gauges. In a small buffer region surrounding the body, 
the form of the perturbation can be found analytically as an expansion for small distances r from 
a representative worldline. Given only a specification of the body's multipole moments, the field 
obtained in the buffer region suffices to find the metric everywhere outside the body via a numerical 
puncture scheme. Following this procedure at first and second order, I calculate the field in the buffer 
region around an arbitrarily structured compact body at sufficiently high order in r to numerically 
implement a second-order puncture scheme, including effects of the body's spin. I also define nth- 
order (local) generalizations of the Detweiler- Whiting singular and regular fields and show that in 
a certain sense, the body can be viewed as a skeleton of multipole moments. 

PACS numbers: 04.20.-q, 04.25.-g, 04.25.Nx, 04.30.Db 
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I. INTRODUCTION 

In general relativity, bodies affect the spacetime geom- 
etry around them, and their motion is, in turn, affected 
by that geometry. Historically, study of this problem of 
motion has focused on the post-Newtonian limit of weak 
fields and slow velocities. In the strong-field regime, bod- 
ies haves typically been approximated as test particles, 
which move on geodesies of a background spacetime that 
is unaffected by them; or, more generally, they have been 
treated as test bodies, which arc accelerated only by the 
coupling of their multipole moments to the background 
curvature. However, gravitational wave astronomy will 
soon allow us to observe binaries with strong fields and 
fast motion, providing fresh impetus to solve the problem 
of motion in the strong-field regime. 

For binaries made up of bodies of comparable mass, 
the problem has been tackled with numerical relativ- 
ity Ql or analytically with effective one body (EOB) the- 
ory [2|. For binaries with extreme mass ratios, in which a 
small mass m emits gravitational waves and spirals into 
a large mass M, the problem has been approached via 
the gravitational self-force approximation: the smaller 
mass is treated as a source of small perturbations of the 
background spacetime of the larger, and that perturba- 
tion forces the small body away from geodesic motion in 
the background Besides modeling these extreme- 

mass-ratio inspirals (EMRIs), the self- force provides an 
essential point of comparison with the methods used for 
intermediate- and equal-mass binaries and can be 

used to fix mass-dependent parameters in EOB P4rij|. 

The gravitational self-force is now well understood at 
linear order in m [H Il2l - fl4| . However, to extract orbital 
parameters from a waveform emitted by an EMRI, one 
requires a theoretical description accurate to second or- 
der, as shown by the following simple scaling argument: 



Suppose z(t) is the body's trajectory and that we have 
an equation for its acceleration a, but that our equation 
differs from the true acceleration by an error 5a. Then 
our result for z(t) will differ from the true position by 
Sz ~ t 2 Sa. Therefore, to obtain a model accurate at 
order 1, we require t 2 5a <C M. An inspiral occurs on 
the timescale over which energy and angular momentum 
change substantially, which is the radiation-reaction time 
t rr ~ M 2 /m, suggesting that we require 6a <C m 2 /M 3 ; 
that is, we need a to be accurate at least through order 
to 2 . This conclusion is also borne out by a more rigorous 
scaling argument [lj| . 

Furthermore, comparisons with numerical simulations 
suggest that the second-order self-force would pro- 
vide a highly accurate description of intermediate-mass- 
ratio binaries and even a reasonably accurate descrip- 
tion of comparable- mass binaries @, Q, both of which 
should soon be observed by Advanced LIGO [H, [l7| . 
The second-order force would also fix EOB parameters 
quadratic in m. Finally, second-order results would allow 
us to determine the error in various first-order approxi- 
mations, such as the geodesic-force approximation used 
in Ref. [ill. 

Hence, there is a need to determine the self-force at 
least through second order, and proceeding to still higher 
order could be useful for studying intermediate- and 
equal-mass binaries. Harte [l9| has shown that even out- 
side the realm of perturbation theory, material bodies of 
sufficiently small size behave as test bodies in a certain 
effective metric, where the effective metric is defined by 
subtracting out the body's 'self-field', which is bound to 
the body and has no direct influence on its bulk motion. 
This split into self-field and effective field overcomes the 
longstanding obstacle in applying Dixon's treatment of 
extended bodies to the self-gravitating case (as discussed 
in Sec. 13 of Ref. [13]). In the context of the self- force, 
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it is a powerful extension of the analogous first-order re- 
sult of Detweiler and Whiting [2l| , who showed that the 
first-order self-force is equivalent (at order m) to geodesic 
motion in a certain effective metric. However, Harte's 
method does not apply to the motion of black holes, be- 
cause it involves integrals over the body's interior. More 
importantly, it does not provide a concrete means of solv- 
ing the Einstein field equation (EFE) to obtain the total 
metric or the effective metric — this must be done in a 
pcrturbative context. 

In recent years, two pcrturbative methods have 
emerged: one that expands the perturbation due to the 
body in a strict power series, representing the body's 
motion using deviation vectors pointing away from a ref- 



erence geodesic |f2|; and one that uses a general ex- 



pansion which self-consistently incorporates the correc- 
tions to the motion, working with an accelerated world- 
line rather than deviation vectors [HI, (3] ■ Both of these 
methods were devised with sufficient scope to be applied 
at any order in perturbation theory. The power series 
method is limited to short timcscales over which the de- 
viation vectors remain small, but there is the prospect 
of combining it with a two-timescale method [I5| to cre- 
ate a valid approximation on the radiation-reaction time. 
The self-consistent approximation can be expected to be 
valid on long timcscales, and it can always be reduced 
to the power series method by performing an additional 
local-in-time expansion. 

Gralla has recently applied the power series method 
to second order [22jj. Shortly before that, I reported the 
second-order results of the self-consistent method [23j | . 
In this paper and its sequel [24|, I present the explicit 
expressions that would not fit within the short letter for- 
mat of Ref. [23} . However, besides recapitulating those 
results, the present paper serves a larger purpose of show- 
ing how the EFE may be solved for an arbitrarily struc- 
tured small body at any order in a self-consistent expan- 
sion. I also show, with slightly less certitude, how the 
nth-order equation of motion will be determined. While 
most of the necessary ideas were sketched in Ref. fl3j . 
here they are put in considerably more concrete form, 
while still avoiding detailed calculations. 

Both the power series and self-consistent expansions 
are based on the method of matched asymptotic expan- 
sions from singular perturbation theory [lj, l25U27j . In 
this method, one uses two expansions, one valid very near 
the body and the other valid everywhere else, and infor- 
mation is fed between them by insisting that they agree in 
a buffer region around the body where they can both be 
expected to be valid. This method avoids modeling the 
body as a point particle (except in a well-defined sense 
discussed later in this paper), since such models are not 
well defined in the full, nonlinear theory. 

A brief review of the method will be worthwhile. Sup- 
pose g M „(e) is an exact solution to the EFE containing 
the small body on a manifold M , where e is an expansion 
parameter that counts powers of the body's mass. Now 
let r be some measure of distance from the body and 1Z 



represent the spacetime's lengthscalcs, excluding those of 
the body itself. For r ~ 1Z, well outside the body, in any 
global coordinates x a in a vacuum region Q, one can use 
the outer expansion 



g^{x a ,e) = g^{x a ) + h^ix"^) 



(1) 



on a manifold Me- (flWi Me) defines an external back- 
ground spacctime with no small body in it, and 



h, lv (x,e)=J2t n h ( $(x;e) 



(2) 



(n) 

describes perturbations due to the body. Each term h^J 
is allowed to depend on e through a dependence on the 
e-dependent worldline 7 that will represent the body's 
motion in Me- This is the defining characteristic of the 
self-consistent approach [l3|, [l4j] ; in the method of Gralla 
and Wald, no such dependence is allowed. As a concrete 
example of the quantities involved in an outer expansion, 
in an EMRI, g^ is the spacetime of the supermassive 
black hole, the external lengthscale is 1Z ~ M, and x a 
might be Boycr-Lindquist coordinates of the supermas- 
sive black hole. 

Now, for r ~ eTZ, very near the body, the metric varies 
rapidly, and there, in coordinates (t,x a ) approximately 
centered on the body, one can use the inner expansion 

g^(t, x a , e) = gx^t, x a ) + Y, z n H$(t, x a ) (3) 

n>l 

on a manifold Mi, where (gi^, Mi) is the body's space- 

(n) 

time were it isolated, and H^' describes perturbations. 
Here the components in the (t,x a ) coordinates are writ- 
ten in terms of rescaled spatial coordinates x a = x a /e. 
These scaled coordinates serve to keep distances fixed 
relative to the body's mass in the limit e<1, sending 
distances much larger than the mass off toward infin- 
ity. The use of a single scaling factor makes the approx- 
imation most appropriate for compact bodies, in which 
the linear dimension is comparable to the mass; I return 
to this issue in Sec. |V1 Scaling only distances, not t, is 
equivalent to assuming the body possesses no fast inter- 
nal dynamics. That is, there is no evolution on the short 
timescale of the body's mass and size. 

In the buffer region around the body, defined by e <C 
r/TZ <C 1, the inner expansion can be expressed in un- 
sealed coordinates and then expanded for r ^> eTZ (i.e., 
for distances that are large on the scale of the inner ex- 
pansion) and the outer expansion can be expanded for 
rCTJ (i.e., for distances that are small on the scale of 
the outer expansion). Since the inner and outer expan- 
sions are assumed to approximate the same metric, the 
results of these re-expansions in the buffer region must 
match order by order in both r and e. References p^.[26j 
contain more details about the underlying assumptions 
in this procedure. But for my purposes here, the essential 
aspect is that an outer expansion can be obtained with 
only minimal knowledge of the inner expansion [I2I [l3| . 
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Furthermore, it can be obtained almost entirely from in- 
formation in the buffer region. Say the nth-order pertur- 
bation is expanded as nju) = J2 P r^ v) in the buffer 
regionQ One must allow negative powers of r, since at 
least part of the field will fall off with distance from the 
body. But there is a lower bound on what the most nega- 
tive power might be at a given order. I allow no negative 
powers of e in the inner expansion, meaning e n h^j must 
have no negative powers of e when written as a function 
of r = r/e, from which it follows that 

e n h$ = ^^'" n) + 0(e n r~ n+1 ). (4) 

Also, since e"/r n is independent of e in f coordinates, it 
must be equal to a term in the zeroth-order part of the 
inner expansion, gi^ u . Say g\ llVl written in terms of the 
unsealed coordinates and expanded, readt@ 

ft" = (5) 

n>0 

We then have hfy = ffw- Therefore, at each order 
in e, the most singular (as a function of r) piece of the 

(n) 

metric perturbation in the buffer region is deter- 
mined by the r 3> m asymptotic behavior of the body's 
unperturbed metric. 

Now, from the perspective of the outer expansion, the 
buffer region surrounds a region B C .Me- If 7 is in S, 
then we can take r to be a radial distance from 7, and the 
form ((3]) will hold true. We then have both the body — in 
the full spacetime — and 7 — in A^e — in the region sur- 
rounded by the buffer, even if B is not diffeomorphic to 
the region surrounded by the buffer in M. (say, if a black 
hole or wormhole resides therein). If mass dipole terms, 
which indicate the position of the center of mass relative 
to the origin of the coordinates, vanish in this coordinate 
system, then the body is appropriately centered "on" 7; 
that is, on their respective manifolds both 7 and the body 
are at the center of the region surrounded by the buffcro 
7 is then a meaningful and accurate representation of the 
body's motion. 

These ideas were discussed at length in Ref. [f3[. In 
the present paper, I focus on their essential practical con- 
sequence: how a solution obtained entirely in the buffer 



1 The coefficients hjl' p are not necessarily independent of r: 
they may have a functional dependence on lnr, as discussed in 
Sec. IIHI However, the terms are still well ordered for small r, 
since r p (lnr) q <C r p (lnr) 9 for p > p' . 

2 The fact that the inner background must be asymptotically fiat, 
containing no positive powers of r, follows from the assumption 
that the outer expansion contains no negative powers of e. 

3 This notion of mass-centercdness based on mass dipole terms is 
not straightforward at high order in e, as discussed in Sec. IHI CI 
I assume that some copacetic notion of centeredness can be im- 
posed. 



region, without reference to details of the body's internal 
structure or global boundary conditions, can be used to 
construct a global solution to the nth-order EFE in the 
external spacetime. By substituting an expansion of the 
form ([Jl into the EFE, one can solve order by order in 
r. Because r is small in the buffer region, spatial deriva- 
tives dominate over temporal ones. Therefore, solving 
the EFE is reduced to a process of solving a sequence of 
flat-space Poisson equations, and the solutions to corre- 
sponding homogeneous equations comprise all the free- 
dom in the general solution. For multipolc number £, 
these homogeneous solutions behave as l/r e+1 and r e , as 
is familiar from classical electromagnetism and Newto- 
nian gravity. The l/r e+1 terms are identified with the 
body's multipole moments or corrections thereto; the r 
terms are identified with free radiation. 

Inspired by this, I define a split of the nth-order solu- 
tion of the EFE into a self-field h^ v and an effective field 
h^ v . Roughly speaking, the self- field is constructed from 
the set of multipole moments and bound to the body. 
The effective field is a vacuum solution to the EFE that 
propagates independently of the body, even though the 
body created it, and in the solution obtained in the buffer 
region, it consists entirely of unknown functions, which 
can be determined only from global boundary conditions. 
My definition of the self-field carries a similar meaning to 
Harte's [l9[ , and the metric g^ v + h^ v defines an effective 
metric similar to Harte's, but because my effective field 
satisfies the vacuum EFE at all orders while Harte's does 
not, the definitions differ at second order and beyond. 

Using the self-field found in the buffer region, one 
can obtain a global solution via a puncture scheme (also 
known as an effective-source scheme), just as at first or- 
der (28l - [33j . Such a scheme, in the present context of a 
small but extended body or black hole, begins by allow- 
ing the expressions for h^ v and h^ v in the buffer region 
to hold in the region B as well. Since doing so does not 
affect the field in the buffer region, it also does not af- 
fect the field values outside the buffer, out in the external 
universe. Continuing the expression for h^ v into B makes 
it into a field that diverges at r = 0: the true self- field in 
the interior of the body, whatever it may be, is replaced 
with this divergent field, and the self-field becomes the 
singular field. But we are not interested in obtaining h^ v 
in the body's interior, nor even in the buffer region, since 
it is already known there. We are interested only in h^ v , 
which is the unknown part of the field in the buffer region. 
Extending the buffer-region expression for h^ v leads to a 
field that is smooth at r — 0, which is the regular field. 
Therefore, in a region covering B we can rewrite the EFE 
as an equation to be solved for h^ v , by subtracting off 
the contribution of the singular field. The singular field, 
or any approximation to it that has the same singularity 
structure, serves as the titular 'puncture' in the scheme. 
At the same time as the field equation is solved for the 
regular field, the 'position' of the puncture is moved via 
the equation of motion for 7. 

With such a scheme, the physical problem in the region 
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covering the body, with all its matter holds, singularities 
(in the case of a black hole), or other oddities (in the case 
of exotic matter), is replaced with an effective problem. 
One well-behaved problem free of singularities (except in 
the case of a black hole) in that region is replaced with 
another, but the variable being solved for is changed from 
the physical field to the effective field h^ v . By appropri- 
ately converting from h^ v to the full field h pv outside that 
region, one can obtain h pv globally in the external space- 
time. One will then have found the physical field outside 
the body solely from a specification of the body's mul- 
tipole moments (which determine the singular/self- field) 
and initial data (which determines the effective/regular 
field), both of which are freely specifiable. 

From this perspective, the various regularization meth- 
ods that have been used in the self-force problem to re- 
move the 'singular part' of the field @ arise only as a 
practical necessity: we cannot determine the physical 
metric inside the body, nor are we interested in doing 
so, which prompts us to replace it with the fiction of a 
singular field solely as a means of calculating the physical 
metric outside the body. 

In the paper proper, I make these notions more pre- 
cise. Section HTl presents the equations to be solved in the 
outer expansion, generalized beyond Refs. [H, [13, HH to 
a wide class of wave gauges. This leads into the anal- 
ysis of the nth-order outer expansion in Sec. IIII1 where 
I show the form of the solution in the buffer region, de- 
fine the singular/self- and regular/effective fields, and de- 
scribe the puncture scheme in greater detail. I also de- 
rive a stress-energy tensor that effectively represents the 
body's composition as a skeleton of (corrected) multipolc 
moments supported on the body's worldline. In addition, 
I show the relationship between the nth-order accelera- 
tion of the body's worldline and the gravitational field 
in its neighbourhood; however, this suffices to determine 
the motion only once one makes an appropriate choice 
of mass-ccntcrcdncss and specifies the evolution of the 
body's multipole moments. 

Section IIVI applies this nth-order analysis at first and 
second order to obtain the first of two explicit results re- 
ported in Ref. (23[: an expression for the second-order 
singular/self- field through order r, which is sufficiently 
high order to find a global second-order solution nu- 
merically via a puncture/effective-source scheme. Pre- 
viously, in Ref. [1 31 ) . the second-order singular field was 
obtained through order r , but explicitly acceleration- 
dependent terms were dropped. The second result re- 
ported in Ref. (23| was that if the body is spherical, then 
7 satisfies the geodesic equation in the effective metric 
+ through second order in e. In the sequel [24| . 
I will provide the detailed expressions involved in the 
derivation of that equation of motion. 

I conclude the present paper in Sec.[V]with a discussion 
of future work and a comparison of my method and re- 
sults to other second-order analyses: that of Gralla [22| . 
mentioned above, earlier results of Rosenthal [HI, HH, and 
a recent discussion of Dctweiler [39 . 



Throughout, I work in units of G = c = 1. Greek 
indices range from to 3. Lowercase Latin indices re- 
fer to spatial coordinates. Uppercase Latin indices de- 
note multi-indices; for example, L = i\---ig. T/l) = 
STF Tl = Tl indicate a tensor symmetric and trace-free 
(STF) with respect to 5 a b- 

II. OUTER EXPANSION IN GENERALIZED 
WAVE GAUGES 

Let D C Me be a vacuum region of the background 
spacctime (i.e., a region in which g^ satisfies the vacuum 
EFE). Furthermore, let it be a globally hyperbolic region 
with Cauchy surface S, such that D = D + (E) and wave 
equations can be solved in D given only initial data on E. 
Now let C A4b be D with the region surrounded by the 
buffer removed; that is, O = D\B. il then corresponds to 
a vacuum region outside the body in the full spacetime. 
This will be the region in which I seek a solution to the 
EFE. 

In ft, the self-consistent method of Ref. [l3[ is mod- 
elled on standard post-Minkowskian expansions in the 
harmonic gauge, as used in post-Newtonian theory [37j . 
It is convenient to use any field variable 4>^ v satisfying 

4>^ = W- v + 0{h 2 ), (6) 

where an ovcrbar indicates trace-reversal with the back- 
ground metric g pv . and the 0(h 2 ) indicates anything 
that becomes of quadratic or higher order when ex- 
panded for small . For example, some candidate 
fields are ^ v = W v , ^ v = W v - ^g^g pa h pa , and 
yiv = H ^v = gltv _ ^pgM". The latter choice is dis- 
cussed in detail in Ref. [38j]; it is the field used in post- 
Newtonian theory in the case that g pv is Minkowski. 

To disentangle the field from the matter degrees of free- 
dom, one may first impose a convenient gauge condition. 
I adopt a generalized wave gauge defined by 

v v 4r = z» vp <t>»p + (7) 

where Z p vp and D M are specified functions of x a and e[3 
I restrict Z^ vp to be smooth on and independent 
of e (in any coordinates that are independent of e), for 
reasons to be made clear in a moment. For simplicity, 
I assume D M to be a smooth function of e, admitting a 
power series D p (x, e) = ^ e n D^(x). I allow D^(x) to 

diverge as l/r" +1 as r — > 0, which preserves the singu- 
larity structure Q of the metric perturbation. 

If cj) 111 ' = h pv , then © is a generalized Lorenz gauge 
condition. If = H plJ , then (JT]) reduces to a general- 
ized harmonic gauge condition when g pu is Minkowski; 



4 One could instead impose SV^" = Z ti vp <t) up + D M , where g V„ 
is the covariant derivative compatible with g^ v , but doing so 
would be equivalent to a redefinition of </>'"' . 
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in the case where <^ v = H' 11 ' and Z^ up = 0, the gauge 
condition can be written as a condition in the full space- 
time rather than in the background [H]. In all cases, the 
gauge condition ([7]) makes the linearized Einstein tensor 
hyperbolic by removing its elliptic part. 

Therefore, with Eq. ([7]) imposed, the exact vacuum 
EFE G^lg] = in f2 is put in relaxed form, as in post- 
Newtonian theory. It becomes a weakly nonlinear wave 
equation, 



(8) 



where 



2R» 



2(Z 



in 



+ g^(Z x prT , x + Z x pa V x ) 



(9) 



becomes the usual wave operator of linearized grav- 



ity when Z^vp — 0, returning E pv f 



hi 1 " 



U<fy 



fUl 



2R» 



p CT 



b pa in that case. The smoothness and e- 
indepcndcnce of Z^ vp on A4e ensures that E^ v pa is a 
nicely behaved operator in the background spacetime. 
The source terms on the right-hand side of the wave equa- 
tion © are 



C u = 2V ( "£>" ) - g^VpDf 



which is a gauge term, and 



rl"' 



i 

8^ 



(10) 



fill 



which can be thought of as the stress-energy tensor of 
the gravitational field, comprising all terms in G^g] 
that are nonlinear in <f> pl/ . Here G^ v [g,(j)] is the full 
Einstein tensor G 1 *" [g + h] after rewriting /i M „ in terms 
of <\> pv , and 8G^[(j>\ is the linearized Einstein tensor 
written in terms of 0**". can be written with the 

gauge condition already imposed or not. Note that un- 
like the Landau-Lifshitz pseudotensor or its generaliza- 
tion to curved backgrounds in Rcf. (38|, contains 
second derivatives of the metric perturbation, meaning 
in addition to its 'gravitational energy' terms, it contains 
terms that should not be so interpreted. 

Equation ([8]) can be solved without constraining the 
body's bulk motion or the evolution of its matter degrees 
of freedom. No stress-energy tensor appears in Eq. ([5]). 
since the small body lies outside f2, which may make this 
lack of constraints seem vacuous. But we can say more 
precisely that in f2, the wave equation can be solved with- 
out constraining the evolution of the multipole moments 
in the buffer region, including that of mass dipole terms 
which shall be used to define the bulk motion. 

Given its smooth dependence on e, C 1 *" can be 
expanded in an ordinary power series C M!y (a;, e) = 
^e"C^(x). I next assume an asymptotic expansion 



(12) 



in which the terms ffiJx; e) are constructed by splitting 
the weakly nonlinear wave equation ([8]) into a sequence 
of wave equations 



Efi u 



\n) 



°(n) 



167rT frT)> 



(13) 



where contains all pieces of t^ v with an overall factor 



ofe" 



o 

(«) 



S 2 G^ [<£(„_!), ^ ( i)]+<5 2 G^[</> (n _ 2) , </> (2 



. . . + 8 n G pv [(j) {1) \, where 8 n G^[<j)] is the nth-order Ein- 
stein tensor. (The notation is described in Appendix El) 
More concretely, at first and second order these equations 
read 



E p.v 



pa 



pa 



LOO _ f-iflV 

V) ~ °(i)' 

^)=^2)-2^[0(i)]. 



(14) 
(15) 



Equations (fTB"]) can be solved in sequence for each <fi^ 
without constraining the motion of the small body, and 
< ^(n)( a: ' e ) ' s taken to be the solution even when the mo- 
tion is faithfully represented (i.e., when the buffer region 
at a distance e <C r/lZ <C 1 from 7 always surrounds 
the body). The dependence on e in each term thus takes 
the form of a functional dependence on 7. I assume that 
despite the dependence on e, each <f>^(x; e) is uniformly 
of order 10 If this is found to be untrue in a particular 
case, then the asymptotic series (fl"2"|) has failed. 

The motion is determined via the gauge condition, 
which ensures that a solution to the wave equa- 
tion (JSJ is also a solution to the EFE (and therefore the 
Bianchi identity). However, I do not wish to split the 
gauge condition |(7J) into a sequence of conditions of the 
type ^ v <t>(n) = Z^vptfn) ~*~^(nV Those conditions would 
imply that • 



v PV(n) 

6^ satisfies the linearized EFE and therefore 
the linearized Bianchi identity, which would constrain the 
acceleration <z M of 7 to vanish [l3[ ■ Yet I do wish to split 
Eq. ||7J) into equations that can be solved for each <ft^y 
To do so, I assume the acceleration can be expanded as 



o"(T,e) = Xy 



a r«)( T ' £ )' 



(16) 



where at time r, each is allowed to depend on the 
past history of the e-dependent worldline 7, in the same 
manner as does. When written in a coordinate 
system centered on 7, all fields involved in the gauge 



condition- 



9p,vi 



Z^ pai and — will naturally de- 



pend on 7's acceleration. So in such a coordinate system, 
one may substitute the expansion for a M , together with 



5 It may also contain a dependence on lne coming from the per- 
turbation of light cones. Such logarithmic terms are discussed 
in Appendix \C\ If they arise, then the statement of uniformity 
should instead read that the coefficient of each (In e) 9 in tfi? < (x; e) 
is uniformly of order 1. 
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the expansion of (f)^ , into the gauge condition to arrive 
at the equation 

n n 

8=1 S=l 

n 
s=l 

for each n, where 5 s returns the sth-order term in the 
expansion with respect to a M and its derivatives. For 
example, S 2 a fJ, a l/ = o^^fo) ~^ a (o)®(2) +^ a (i)^('i)- " m di- 
cates evaluation at a M = a^ y Equation (JTTJ) can also be 
written as an equation for <f>^ in terms of the lower-order 
fields: 

There is no obvious way to express the 6 s operation 
in a covariant manner, and one could do away with it 
entirely by instead solving approximate, rather than ex- 
act, equations for each <pf^y In 'little o' notation, these 
equations would read 

e(v^ - ZV<f) - £>?!)) = o(e), (19) 
E e " ( V »C) - - D (n)) = °( f2 )' ( 2 °) 

n=l,2 ^ ' 

etc., where o(e") means terms that approach zero faster 
than e™ for small e. Equations and (f2"U|) follow di- 
rectly from Eq. ([7]) , given the assumed uniformity of each 
4>(n\ ■ However, the expansion of the acceleration helps 
to systematize the calculation, and it trivially leads to 
order-reduced equations in cases where derivatives of the 
acceleration would appear on the right-hand side of the 
equation of motion. 

As in traditional post-Minkowskian expansions, these 
gauge conditions are closely related to the conservation 
of the source T pv . Taking the divergence of the wave 
equation (J8j) and using the identities \7 U CP 1 ' = U\D P and 
V„E^ pa <f>P° = □ (V„</>^ - Z» vp r p ), we find 

Off" = -167tV„t^, (21) 

where I have defined 

= V v <ff v - Z^ vp( j) vp - D^, (22) 

in terms of which the gauge condition reads = 0. 
From this it follows that if (i) V u t^ = in tt and (ii) 
= V n H^ = on the boundary <9f2, where n a is a vec- 
tor normal to dQ, then H p = everywhere in £1. Equa- 
tion together with the Bianchi identity shows that 



if G^[g,(j)] = o(e n ), then V„t^ = o(e n+1 ). Combining 
these results, we can show by induction over n that if 
(i) the data on the part of £ outside the buffer region 
satisfies the gauge condition, (ii) the gauge condition is 
satisfied in the buffer region to an accuracy H p = o(e n ), 
and (iii) the first n wave equations are satisfied, then the 
gauge condition is satisfied up to o(e ra ) errors everywhere 
in Q (putting aside issues of secularly growing errors in 
H 11 ). Usually, one has that the conservation equation 
Vi/(T Ml/ + t^ v ) = together with gauge-satisfying ini- 
tial data suffices to enforce the gauge condition every- 
where; here, since we exclude the region containing mat- 
ter, enforcement of the gauge condition in the buffer re- 
gion takes over the role of conservation of material stress- 
energy. 

Note that in previous work [lU, [l4|, [Hj], I used the 
metric perturbation h pv as the field variable, rather than 
something based on its trace-reverse. Had I used a field 
§p,v = hf^v + 0(h 2 ) as the variable in this section, ev- 
erything would have carried through almost identically. 
The gauge condition would have involved the divergence 
of the trace-reversed field, {g*g ua ~ ha^^^^Xo = 
Z tI vp (j) U p + Dp, and the sequence of wave equations for 
(n) 

<ppa would have followed from an expansion of the exact 
EFE i? M „[g + h] = 0, since the linearized Ricci tensor in 
a wave gauge takes the form of a wave operator acting 
on in the same way the linearized Einstein tensor 
takes the form of a wave operator acting on h^. In the 
present paper, I favor the trace-reversed perturbation for 
the sake of the simpler gauge condition. 

Also note that throughout this section, I have used 
the language of post-Newtonian theory, in which the ex- 
act EFE is said to be relaxed by adopting a gauge in 
which the EFE can be solved without constraining the 
source. One then seeks a field 4>p v that solves both the 
relaxed EFE and the gauge condition as accurately as 
possible. In self- force literature [3, [ID, [22[, there is a 
notion of relaxing the gauge condition, meaning one im- 
poses V„0^ = but then allows V„e</>^ = 0(e 2 ), which 
avoids constraining the motion at order e. These two no- 
tions of relaxation are similar but distinct, and there is 
no sense in which I relax the gauge condition on a pertur- 
bation: \7 u (j)^ = is simply never imposed; only within 
the context of a Taylor series (i.e., regular perturbation 
theory) would the condition V„<^" = imply the condi- 
tion V„0^ = 0, which would then have to be relaxed in 
order to self-consistcntly evolve the position of the body 
and the field it produces. 



III. nTH-ORDER OUTER EXPANSION 

Before constructing explicit solutions to the field equa- 
tions for a particular choice of , Z% CT , and D M , I 
discuss the general form of the solution at any order. 
This involves two steps: first, deriving a local expres- 
sion for a certain self-field in the buffer region; second, 
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using that local expression to obtain the global solution 
via a puncture/effective-source scheme. I first derive the 
general solution to the nth-order wave equation (|13p in 
the buffer region, subject only to the matching condition, 
with neither the body's composition nor other boundary 
conditions specified. 

From the form of the field, I also derive an effective 
form for the body's nth-order stress-energy tensor in 
terms of the body's multipole moments and the correc- 
tions to them. Finally, from the gauge condition (fl7|) in 
the buffer region, I derive the relationship between the 
nth-order acceleration of 7 and the metric perturbation 
in its neighbourhood. 



A. General solution in buffer region 

In the buffer region, I work in local coordinates (t, x a ) 
centered on 7 and use 6 a b to raise and lower lowercase 
Latin indices. Any choice of local coordinates will do, so 
long as the metric reduces to Minkowski on 7; that is, 
9ftv = Wftv + 0{r), where r = y / %Jx^a? is the spatial dis- 
tance from 7. For the explicit results in Sec. IIV1 I adopt 
Fermi- Walker coordinates, in which t is proper time on 
7 and r is the proper distance along a geodesic perpen- 
dicular to 7. But prior to that, excepting a brief detour 
in Sec. IIIICl I leave the coordinates unspecified. 

I write 4>^n) as ^(n)^! r i ^)i where the unit vector n l = 
x % jr. Favoring n a over the usual angles (8, cf>) is useful 
because the smooth background metric g^ v will naturally 
be a power series in x a = rn a . I assume that 4^,) r ' 
possesses an expansion of the form 



: 9m a 

E 

p>-n 1=0 <}=0 



EE^ h -)%;V^> ( 23 ) 



where the coefficients 



time, and Ul 



9 {npq£) 



are smooth functions of 



ni t \ is STF with respect to 5 a b 



The decomposition in terms of h j is equivalent to an 
expansion in spherical harmonics [39|]; standard useful 
identities involving this decomposition are provided in 
Appendix [B] The lower limit on the sum over p, which 
makes r~ n the most singular power of r, comes from the 
matching condition, as noted in Sec. |U The logarithmic 
terms arise from the perturbation of light cones. One 
can expect the solution to the exact EFE to propagate 
on (and within) null cones of the exact spacctime, and 
given that the mass of the body induces a logarithmic 
correction to the retarded time, logarithmic corrections 
then naturally appear in . This effect is well known 
from solutions to the EFE in harmonic coordinates (see, 
e.g., Rcfs. [H,|43l). Other causes may also lead to loga- 
rithmic terms. For generality, I allow logarithms at any 
value of n, but I assume that for each finite n, p, and £, 
the highest power of lnr is a finite number q m n K (n,p,£). 
It will be seen that this maximum is required to obtain a 
solution. For simplicity, to make sure that term-by-term 



differentiation is valid without worrying about issues of 
convergence, I also assume for a given, finite n and p, £ 
has a maximum ^max(^iP)- 

I proceed by substituting this expansion into Eq. (|13j) 
and solving for the coefficients 4>(npq£)(t)- To facilitate 
that procedure, I split the wave operator into two pieces, 
E^ pa = gjjgZd'di + AE pv pa . The first operator is a 
flat-space Laplacian, which reduces a term of order r p 



to a term 
operator, 



r p-2 



(treating r p (\nr) q ~ r p ). The second 



AE^ pa = e ^p° - 9^9% 



(24) 



takes a term of order r p and returns terms of order r p 1 
and higher. I break this down further into 

A£^V= E (P ' )A ^V. ( 25 ) 
p'>-i 

where ( p ^AE^ V pa takes a term of order r p and returns a 
term of order r p+p . The wave equation (|13|) can then be 
written in the form of a Poisson equation, 



d l d % 



_ pp,v 



167^ 

(«) 



(26) 
(27) 



Because t?\ is made up of nonlinear combinations of 
^{n') ( wnere n' < n), it necessarily has an expansion of 
the form 



E E^-^Vg^. (28) 

p~> — n q.£ 



Similarly, because Aff" prr4> p (Z\ is simply a combination 
of unit vectors and derivatives acting upon <j> P fZ\, it has 
an analogous expansion. Assuming can likewise be 
expanded as 



°(n) 



p> — n q,£ 



(n,p,q,£) 



(29) 



the source in the Poisson equation then also has an ex- 
pansion 

p h = E EE ^^v^La** ^ 

p>~n £=0 9=0 

where C = C(n,p) and Q — Q(n,p, £) are finite for finite 
n, and the coefficients are given by 



P, 



fj,i>L 



(lnr)"" Yl 



E 

p'—~n I' ,q' 



°(»,p, 9 ,£) 



(n,p,q,£) 



L 
(Q) 

(31) 
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Here [/]^ denotes the coefficient of (lnr) 9 ?!^ in the series 
/ = J^i' <7 / (h ir ) 9 [/] i ^L't where the coefficient [f] L is 



given by Eq. (|B2|). 

Because is an eigenvector of the flat-space Lapla- 
cian, satisfying r 2 d l difiL = —£(£ + l)nz, (see Ap- 
pendix [Bj, one can easily evaluate the left-hand side of 
Eq. (f2"7) . Doing so, and substituting the expansion of 
P/l^ into the right-hand side, leads to 



(n) 



x{b(P+l)~^ + l)]C^) 
+ ( q + l)(2p+l)^ L pq+ie) 

+(?+l)(« + 2)C£«+M>} 

= E EE^(in,)M: 

p>-n e=o q=0 



L 

(npqi) ' 



(32) 



where I have defined 



This 



^npqi) = for q > q max . 
can now be solved by equating coefficients. Since the 
coefficient of r p ~ 2 in P?\ is constructed from 6^,-, for 

(n) V(n') 



< n and 



ipvL 
(n,p',q,i) 



for p' < p, & term 



(n,p,q,l) 



Oil 



the left-hand side is sourced by lower-order terms on the 
right-hand side, allowing one to solve order by order. 

The form of the solution obviously differs depending 
on whether p(p + 1) — £(£ + 1) vanishes. I define 



-p - 1 if p < 0, 
p if p > 0, 



(33) 



which is the unique non-negative solution to p(p + 
1) - £ p (£p + 1) = 0. I call 0^0^) a homoge- 
neous mode, inspired by the fact that rPn L p 4*^pot ) 
satisfies the homogeneous Poisson equation. This can 
easily be seen by ignoring the logarithm terms, in 
which case the homogeneous Poisson equation reads 



r p-2 



n L [p(p+l)-£(£+l)}^ M) 



0. 



Now, for a given p, first consider the modes for which 



^ £ p . Each mode 



modes 



,pvh 



(npql) 



and 



is coupled to the two higher-g 

,J {n,p,q+l,i) ^(n L p,q+2J)- Start Wlth 1 = ^max, 

in which case the two higher-g modes vanish by definition; 
this is the reason why one must assume finite g max in 
order to obtain a solution. If g max > Q, the Poisson 
equation reads [p(p + 1) - £(£ + 1)] 4^qi) = °' implying 

tfnpq^l) = - !t fo H° WS tliat for ^ &p, 9max = 

and the q = q ma ,x term in the Poisson equation reads 

[ P (p+i)-£(£+i)]^ Qe) =pi: p L Qey 



or 



ipvL _ 
(n,p,Q,i) ~ p(^p 



Pi 



pvh 



(npQi) 



1) 



1)' 



(34) 



For q = Q — 1, one then immediately finds 



ppvh 

( n iPiC — i,' 



(n,p,Q-i,e) - p( p + I) - £(£ + 1) 

Q(2p + l)PZ L c 



(npQl) 



[p(p+l) -£(£+!)}' 



(35) 



These two results are sufficient to find <p^p q ^ for any 
< q < Q - 2, since 



"(npql) 



1 



p(p+l)-^ + l) 

-(0+1)(2j>+1^" £ 



(npqi) 



(n,p,q-\-l,£ 



l)(? + 2)C, +2 i) 



(36) 



Now consider modes for which £ = £ r 



mode 



nivL 
{npql) 



is coupled to one higher-g mode, 



In this case each 

itvL 



(n,p,q+l,. 



Following the same procedure as for £ ^ £ p , one finds that 



for 



is at most Q + 1. For q = Q + 1, 



(",P,S+1,- 



Pi 



pvL p 



(npQ£ p ) 



(e + i)(2p+i) - 



(37) 



Note that this raises the power of the highest-order loga- 
rithm to one higher than that in the source. If the source 
contains no logarithms, logarithms arise in the solution 



only when P { 
l<q<Q, 



pvL v 
(npOlp) 



7^ 0. Next, for any q in the range 



iP-vLp 
^(npqip) 

Last, for q = 0, 



^pvLp 

(n,p,q-l,e p ) _ q + 1 iUvL p 

(2p+i) 2p+v ( - n 'P^ +i ^y 



q 



OLvLp 

(n P oe p 



is undetermined. 



(38) 



(39) 



When the general solution is constructed in this way, 
in it is either a homogeneous mode 



each coefficient 



(npql) 



or directly proportional to a mode Pu^f* in the source 



P, 



(«)• 



Since P, 



modes 



,pvL 



pLlsL 

npqi) 



is constructed from the lower-order 



(n'p'q'i') 

that each ^ vL 



with either n' < n or p' < p, it follows 
(npqi) ^ S e ither a homogeneous mode or con- 



structed from combinations of homogeneous modes (or 
constructed from a mode of C^ v ). Therefore, the set of 



ipvLp 
(npSlp) 



comprise all the freedom in 



homogeneous modes 
the general solution. 

The homogeneous modes can be divided into two cat- 
egories: those for p < and those for p > 0. For p > 0, 
they behave as r , and they (combined with the higher- 
order-in-r terms constructed directly from them) corre- 
spond to free wave solutions to the homogeneous wave 
equation, not directly influenced by the presence of the 
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body. If C^ u vanished and the buffer region did not sur- 
round a body, these modes would fully characterize the 
perturbation. 

For p < 0, the homogeneous modes are the coefficients 
in a standard large-r multipole expansion far from a 
source, as arise in electromagnetism and post-Minkowski 
theory [3^, |4l| . A 1/r homogeneous mode is a monopole; 
a 1/r 2 homogeneous mode, a dipole; and so on. I there- 
fore define the (unnormalizcd) multipole moments 



D (n,—£-l,(M) 5 



(40) 



solving order by order in r and lnr in sequence (start- 
ing at the most singular order in each case). III. The 
highest power of lnr in the solution is increased when 
the source in the Poisson equation (f2"T|) contains a mode 
P{npqi )' wnere i s given by Eq. ((55)) . IV. Each term 
in the general solution is determined either by the source 
C? , — 167tt^ , by the body's multipole moments or grav- 
itational corrections to them, or by a homogeneous mode 



which are symmetric in fj, and v and STF in i\ . . . ig. 
In the solution (|23[) . these moments appear in terms 
■^(n) ^l/t^ +1 . Each moment (except the monopoles), can 
be split into a mass and current moment , 



'(») 



M, 



(n) 



(41) 



where the mass moment M?K L is the even-parity part of 



ml 



(n) 



t(n) , satisfying 



M, 



ft(jii)i2—H 
(n) 



moment is the odd-parity part, satisfying S^' 



S 



fi\jii]i 2 —ie 
(n) 



and the current 

(n) 



Since the most singular term in h$, hfy "'/r™, is 
equal to the term g[™l/r n in the large- r expansion of 
the body's unperturbed metric, the multipole moments 



I(n) N 1 that appear in tif^j 



are related to the mul- 



tipole moments of the inner background spacetime (de- 
fined by, e.g., Geroch [iU and Hansen [43|)i£| Therefore, 
specifying the moments 1^ ~ x is equivalent to specify- 
ing those of the body. The remaining moments, I^\ L 
for £ < N — 1, are gravitationally induced corrections to 
those of the unperturbed body. They will be discussed 
further in the following section. 



,fiuL v 



for 



This means that the homogeneous modes r( np(W 
p < are determined by the presence of the body. More 
precisely, the moments I^-^ -1 are fully determined by 
the matching condition, which plays the role of a bound- 
ary condition. The other moments, for £ < n — 1, 
will be found to involve these fully determined moments. 
Contrariwise, the homogeneous modes for p > are en- 
tirely undetermined by the matching condition; they are 
determined only once further boundary conditions are 
imposed. 

We now have the following conclusions: I. In the buffer 
region, a solution of the form (|2"3")l to the nth-order wave 
equation (|13[) can be obtained to any order in e and r. II. 
The solution may be constructed by a simple algorithm, 



B. Minimal solutions, singular /self- fields, and 
regular/effective fields 

Given that the homogeneous modes <jtf vLp n „ ^ contain 
all the freedom in the solution to the wave equation, the 
solution can be naturally split into several pieces. 

Consider a wave equation E^ v P o4> P ^ = F^ v , with some 
source F^ u . I define the minimal inhomogeneous solution 
<^(n) PCT to this equation to be the particular solution in 

which all homogeneous modes vanish, (ff^^o i ) = ®- ^ 
contains the minimal number of homogeneous modes of 
any particular solution. 

Now take some other particular solution <fij£) ^° ^ e 
same equation. It differs from 4> l ^ a in that it con- 
tains nonvanishing homogeneous modes, as determined 
by some choice of boundary conditions. Consider the 
homogeneous equation E^ pai'fW) = 0- I define the min- 
imal pth solution to be the solution containing 

^(n^o i ) as ^ S oruv nonvanishing homogeneous mode, 

livL , 



with 



WAV) 



for all p' ^ p. We then have 



(n) 



(«) 



E 



(p) 



(42) 



where in each 



(P)^(n)' ^(n^.O^p) 



takes the value it has in 



4\n) ■ ^ n ^ nc simple case where no logarithms ariscQ one 

has (p^^for.ri) = T, p >> p rP ' { P )<t>(ny) (*>")> where the 
lowest-order term is the homogeneous mode itself, 



(43) 



The Geroch-Hansen moments are defined only for stationary 
spacctimes, while gi M „ may depend on time for a physically re- 
alistic system. However, the time dependence is slow compared 
to the body's lengthscales. So the definition may still be used. 



7 I have proved by brute-force analysis of E^ v pcr that any minimal 
pth solution contains no logarithms if p > — 1, though I omit the 
proof due to its tcdiousncss. I hypothesize that the same holds 
true if p < — 1, but I have not proved that to be the case. 
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and the higher-order terms are 



rP '( P )C,P') 



= EE 

£ p"—p 



-1 \(p'-p")J^^ v 



if"? 



pa' (p) p ") 



l> L 



p(p + 1) - £(e + 1) 



(44) 



where [f] L is given by Eq. (|B2j) . 

With these preliminary definitions, I now define an 
nth-order singular/self-field (j>^ and regular/effective 



field 



lR.uv 



, along with their sums 



\n) > 
»(„) . 



(45) 
(46) 



as certain pieces of the physical solution to Eq. (113)) . By 
'physical solution' I simply mean (|23p with all the coeffi- 
cients determined. In other words, given a metric in the 
buffer region around the body, extract the field <^ v from 
it and decompose that field into the form (f2"3"|) ; I define 
^Suv an( j ^Ra^ ag certain pieces of the result. 

Every term in the nth-order solution to Eq. (|13[) is a 
combination of homogeneous modes (and modes of C^). 
I wish the regular field to contain everything in the so- 
lution (constructed in the manner of Sec. IIII A|) not di- 
rectly involving any of the body's multipole moments. 
At first order, where is the only source in Eq. ([TBI . 

this leads to the simple definition that (j)^" is the sum 
of minimal pth solutions for p > 0, with the homoge- 
neous modes therein taking the same values as in the 
physical field 4>^"y Note that this </>^g" is a smooth 
field at r — 0. At the first two orders in r, it reads 

<t>{i) = ^(iao.o) + r ^a?i,o,i) n « + °( r2 )- Beyond first or- 
der in e, one could define the regular field in the same way, 
as the sum of nth-order minimal pth solutions for p > 0. 
Again, this would be a smooth field at r = 0. However, 
in the spirit of Detweiler and Whiting [2ll |. I wish the 
total regular field (j} R P" to be a solution to the homoge- 
neous EFE. To that end, I define dfi't 1 ' to be the sum of 

; (n) 

(i) all nth-order minimal pth solutions with p > (again, 
with the homogeneous modes taking the same value as 
in the full field ) and (ii) the minimal inhomogeneous 
solution to 



where t^K [4>] is as defined in Sec. |TTJ It follows that 



(47) 



iR/ii/ 



satisfies 



i>(j 



= —VOTXT 



(48) 

Since every function involved is smooth at r = 0, <p R P L ' is 
a smooth solution to the homogeneous wave equation not 



only in the buffer region but also in the region BcMe 
that is surrounded by the buffer. Note, however, that 
my definition of <f R L il/ is local, in terms of a certain series 
solution in a certain coordinate system. 

I next define the nth-order singular/self-field to com- 
prise everything else in 4>^n)'i ^ n ^ s means that in the so- 
lution to the nth-order wave equation constructed in the 
manner of Scc lIII Al the singular field contains every term 
directly involving at least one of the body's multipole 



moments. More precisely, 



(n) 



consists of the sum of 



all nth-order minimal pth solutions with p < 0, together 
with the minimal inhomogeneous solution to 



,Sp<T 



°(«) 



16tt(t^ 



(n) '(n) 



(49) 



Again, the homogeneous modes defining the minimal pth 
solutions are given the values they have in the full field 
4ff^y From that full solution, the singular/self-field con- 
tains all the terms that diverge at r = (along with the 
higher-order terms that follow from its definition). At 
first order, where t^X — = t^[^> r ], <fifi) i s a solution 
to the governing wave equation (| 13[) in f2. At higher or- 
ders, it ceases to be, and it begins to directly involve the 
lower-order regular fields. 

With these definitions, I have split the nth-order field 
in the buffer region into a singular and a regular piece, 



(n) 



\n) > 



and, likewise, the full solution has been split as 
4F V = 4> s ^ + <t> K ^. 



(50) 



(51) 



As desired, cj) K P L ' shares the properties of the Detweiler- 
Whiting regular field. It is a smooth solution to the 
homogeneous wave equation even at r = 0. As a ho- 
mogeneous solution, it propagates independently of the 
body and is determined by global boundary conditions 
that specify, for example, the incoming wave content in 
the full solution. The singular /self-field is the rest of the 
full field; that is, it is the part that docs not propagate 
independently, instead being bound to the body. 

All of this puts aside enforcing the gauge condition. 
Ideally, one can impose separate conditions on the singu- 
lar and regular fields: 



DP, 



(52) 
(53) 



This 
ZP, 



ensures the total field 4>p u satisfies = 
P a4> p(J + DP, but it puts all the dependence on DP 
into the singular field, in accordance with how C pv ap- 
pears only in the wave equation for the singular field. In 



Eqs. flU} and gSJ, t ( ^[0 r ] is written with the gauge 

condition on <p R P 1 ' already imposed, which ensures that 
(j) R P L ' is a solution to the full vacuum EFE, not merely 
its relaxed version. 
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However, there is a potential difficulty in imposing 
gauge conditions separately on the singular and regu- 
lar fields: since the singular field contains none of the 
minimal pth solutions with p > 0, it may not have suf- 
ficient freedom to satisfy its gauge condition. If this is- 
sue arises, one may abandon the goal of imposing the 
separate gauge conditions ([5^1) and (|53[) on each of the 
fields, thereby abandoning the goal of having (j^v sat- 
isfy the full vacuum EFE. Alternatively, one may incor- 
porate some particular minimal pth solution (or however 
many are necessary) into </> S/il/ that ensures it satisfies its 



gauge condition. 



,Rpi/ 



will then contain the minimal pth 



solutions with <?^p P p ) replaced by 4>^^ p 



whcr 



(n,p,0,p) 



p) 



(n,p,0,p) ' 



is the homogeneous mode introduced into 



(j>^) in order to satisfy the gauge condition. ^V^, JV) 
will not be uniquely specified in that case, since, among 
other reasons, one could always add to it another mini- 
mal pth solution of the same p but contributing nothing 
to Eq. ([52)1 . Therefore, this will make the split into sin- 
gular and regular fields non-unique. But any convenient 
choice may be made, and </> RA " y will remain a smooth so- 
lution to the vacuum EFE in B and in the buffer region 
surrounding it. In any case, this issue does not arise at 
the orders explicitly considered in Sec. IIVI 



C. Gauge condition, corrected multipole moments, 
and equation of motion 

This consideration of the gauge condition leads well 
into a more detailed discussion of satisfying that condi- 
tion. As shown above, the homogeneous modes contain 
all the freedom in the solution of the form (f2"3")l . And from 
the perspective of the wave equation (fT3"]) . these homoge- 
nous modes can be arbitrary (smooth) functions of time. 
However, they are constrained by the gauge condition. 

We can analyze this condition in a manner similar to 
the analysis of the wave equation. First, note that 



V,<) 



d b<t>(n) 



Or 



(n) 



V vp 



(n) 



(54) 



After substituting the ansatz (|23l) and making use of the 
identities (|B14[) and (|B15[) . the first term becomes 



d b<t>t) 



J2 rP-^lnry 

p> — n 



+(i-M(?+i)<i! g+M „ 1) 

^ (£+l)(£ + p + 2) ijL 

0j i^(n,p,q^+l) 



2^ + 3 
21 + 3 



n L . 



(55) 



Since both Z^ vp and the Christoffel symbols of the 
background are smooth, they can be expanded in Tay- 



lor series Z M 



vp 



E^o^VsU and 



S s >o r Tn r tp,s | r =o' Puttin S tliese results together, one 
finds that the nth-order gauge condition (fTTj) becomes 



= £ 5'\ (1 - S m )( P - i + l)^,,,^) 



s=0 



+ (l-5,o)(g + l)< z : > s . p . g+1 . £ _ 1) 
^ (l + W + P + 2) ijL 



21 + 3 



l 3^(n-s,p,q,£+l) 



(l + l)(g+l) ljL 

2£ + 3 % ^(n-a,p,q+l,t+l) 

p-1 /pp I p(T „p _ <7U \ 



p' — — n V 



x n L 



P-P'-livpL' 



(p-p'-l)! 



(n— s,p',g,^') 



jiL 



^t$(n— s,p— l,t/,^) ^ (n— s,p,q,^) ( ' 



(56) 



where [f] L indicates the coefficient of fiL in the expansion 
/ = J2if] L ™L, given by Eq. (|B2|) . Here we have an equa- 
tion for tfZ* tqA and ^ L p + q 2 e+2) in terms of <^^ q+lA , 
^(n L pq\i 1+2)1 and terms of lower n and/or p. 

Equation (|56p is solved in the same manner as was the 
Poisson equation (|32p . proceeding order by order in r, 
beginning with the most singular terms, <P^_ n q ^ ■ Since 
the gauge condition couples modes of degree I to those 
of degree £ ± 2, doing so is not quite as straightforward 
as in the case of the wave equation. However, it may still 
be readily accomplished. 

Because the homogeneous modes contain all the free- 
dom in the solution, I am interested only in how Eq. (|56[) 
constrains those modes. It could be the case that the ho- 
mogeneous modes provide insufficient freedom to satisfy 
the gauge condition. However, one may always impose 
any desired generalized wave gauge, since the transforma- 
tion from any other gauge can be constructed by solving a 
sequence of wave equations for the gauge vectors^ Since 
Eq. (|23[) is a very general form, it is very unlikely that it 
would be unable to satisfy the gauge condition, most of 
which is satisfied automatically by virtue of the Bianchi 
identity. 

So let us examine how Eq. (|56p constrains a given 
homogeneous mode. The equation yields at least three 
constraints on each homogeneous mode: For q = and 



In Ref. Il4l . I showed that the Lorenz gauge condition could be 
imposed at all orders by solving a sequence of wave equations. 
There, I showed that if the condition could be imposed in the 
buffer region, then it could be imposed everywhere in the external 
universe outside the buffer region. It follows from the results of 
Sec. IHlXl that the condition can be imposed to arbitrarily high 
order in the buffer region by solving the relevant wave equations, 
which completes the proof. 
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iu(iL p ) 
J (n,p,0,e 



The 



I — 1 = £ p) the first term in Eq. (|56|) involves < 

a certain part of the homogeneous mode <^p P ( y 
equation is then one for this quantity in terms of lower- 
order modes and three modes of the same n and p, 
aivL p +2 < h^ vL " anc | dtf vL "? r J . For q = and 

^(n,p,0,« p +2) ' ^(n,p,l,£ p )' "(n,p,l,£ p +2) 1 

£ + 1 = £ p , the third term is a different piece of the same 
homogeneous mode, L p F g j , and we arrive at an equa- 
tion for that quantity. For q = and £ — £ p -\, we have 
an equation for the time derivative of the homogeneous 
mode. Therefore, the tensorial form of the modes will 
be constrained by the first two equations, and its time 
dependence by the third. 

The constraints on the homogeneous modes appearing 
in the regular /effective field are not of particular interest, 
since the regular field will be calculated numerically from 
a puncture scheme. The constraints on the homogeneous 
modes appearing in the singular /self-field, however, are 
essential, because the singular field will be used to define 
the puncture, the form of which must be determined an- 
alytically. Furthermore, the homogeneous modes in the 
singular field are the multipole moments ■ They con- 
tain information about the body and its interaction with 
the geometry, and in the case of a material body, the con- 
straints on them should correspond to the conservation 
of the body's stress-energy tensor. 

The multipole moments come in two types: the body's 
unperturbed moments, which behave as if the body were 
isolated, and corrections to them. As mentioned in 
Sec. IIII Al the moment occurring in the most divergent 
piece of <j)^ at a given n is an unperturbed moment, 
corresponding to the I = n — 1th moment of the inner 
background spacetime g\p V . This moment may be freely 
specified modulo the constraints placed on it by the gauge 
condition. For example, in the case Z^ vp = = 0, 
Eq. ([56]) with q = 0, p = — 1, and £ — 1 constrains the 
first-order monopole moment to have the form (l3j 



(57) 



where the matching condition implies m is the ADM mass 
of gi^y ; and Eq. ([56]) with q = 0, p = 0, and £ = con- 
strains m to be constant, but its value is freely specifiable. 

For £ < n — 1, I?K L is a correction to the body's £th 
moment. From the perspective of the wave equation, 
this could be freely set to zero, effectively incorporat- 
ing the correction into the body's unperturbed moment. 
Therefore, the finite corrections come entirely from the 
gauge condition. A consequence of this is that the cor- 
rections, as obtained entirely within the buffer region, 
without specifying global boundary conditions, will al- 
ways directly involve the unperturbed moments, as men- 
tioned in Sec. IIII Al terms in the gauge condition involv- 
ing only the homogeneous modes with p > cannot have 
negative powers of r, and therefore cannot constrain the 
modes with p < 0. The corrections to the multipole mo- 
ments may arise from the nonlinear effect of couplings 
between multipole moments or between a multipole mo- 



ment and the regular field. They may also arise as the 
linear effect of coupling between a multipole moment and 
the external curvature. However, these corrections will 
not be uniquely determined, because one can always add 
a minimal pth solution satisfying the homogeneous equa- 
tion 5°(Vv(f>^ — Z^^p(j) u ^) — 0. That is, one can always 
add a field indistinguishable from that due to the body's 
unperturbed moments. Such additional terms may be 
freely specified. 

In addition to these constraints on the multipole mo- 
ments, the gauge condition also determines an equation 
of motion for the body, in the sense of providing a re- 
lationship between the acceleration of 7 and the field 
in the buffer region. For simplicity, consider the case 
Z^vp = = 0. Also, take the background metric to 
be 

9nv = f)iiv — c laiX l u il u u + 0(r 2 ), which is the form 
it takes in Fermi- Walker coordinates. Now, with these 
simplifications, consider the contribution of cf>^ to the 

gauge condition. Ignoring the 5, the £ = 0, p = — 1 piece 
of Eq. ([SB]) yields 



1 

-in 



4d t mu a + ^ma a (58) 
Au f3 V fj {mu a ), (59) 



where the integral over a sphere picks off the £ = piece, 
and I have made use of Eq. ([57)) . Equation ([59]) tells us 
that to the p = — 1, £ = part of the gauge condition, 



the first-order field 



/,/"' 
(i) 



contributes what looks like the 



left-hand side of an equation of motion. The right-hand 
side of that equation of motion will be contributed by 
higher-order fields. Write the nth-order gauge condition 
(fI7]l as 



(60) 



Picking out the p = — 1 term using Eq. (|56p. integrating 
over a sphere to pick off the £ — piece, and making use 
of Eqs. {5SJ and ([B18]l . one finds 



4ma 



(n) = _ E ^" " \ o^'( 2 ^K + l,0,0,l) + ^ '(n> '+1,0,1.1)) 



+ ^&+l ) +« MJ K+l)+^P+l) 

-2 L-^J 

+ E E 



(-J/-1)! 



<2p")l(-2p' - 2p" - 1)!! 
p<=-0'+i) P "=o y ' v y F ' 

xa il ii ... v ,v' ( _ P ,_ 1 _ 2 p» ) (r^+r^) 



x aM-p'-^p") 

A ^(n'+l,p',0,-p'-l-2p") 



(61) 



where I have used the first-order gauge condition to set 
dtm = 0. The full equation of motion then reads ma^ = 

m (af ) + ea (i) + f2a f2) + •••)• 
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To turn this relationship into a useful equation of mo- 
tion, one must actually solve the wave equation and 
gauge condition to sufficiently high order to determine 
the content of the coefficients <fi^ L p q f] appearing in 



Eq. (|61[) . Furthermore, one must impose some condition 
on the worldline in order to ensure that it accurately rep- 
resents the body's motion. This will take the form of a 
condition on the mass dipole moments M^ 1 . 

The first mass dipole moment that appears in </> A " y ', 
M^ 1 , corresponds to the mass dipole of the inner back- 
ground spacctime. It is a measure of the 'position' of the 
unperturbed mass m relative to the origin r = 0; a trans- 
lation x a — > x a + Sx a (t) applied to the m/r monopole 
induces a M^ni/r 2 dipole. (Less intuitively, a transfor- 
mation x a — >• x a + 8x b (t)n a b also induces a mass dipole.) 
Equivalently, the mass dipole may be shifted by adopting 
a slightly different worldline about which to center the co- 
ordinates. Hence, this mass dipole is set to zero to ensure 
that 7 represents the body's bulk motion. Analogously, 
the corrections to the mass dipole, M? 1 ? for n > 2, are 
measures of the 'position' of the corrections to the mass 
m. In other words, in the local coordinates (t,x a ) cen- 
tered on 7, the mass dipole terms are a measure of the po- 
sition of the center of energy, including 'gravitational en- 
ergy', relative to 7. The lowest-order mass dipole seems 
to be on firmer ground as a measure of center of mass 
than are its corrections as measures of center of energy, 
since it can be related to a well-defined notion of the 
mass dipole of the inner background spacetime. But for 
an observer sitting in the buffer region, the corrections 
to the mass dipole are indistinguishable from the unper- 
turbed mass dipole, since they take the identical form 
in the metric. More strictly speaking, a certain piece of 
the corrections is indistinguishable from the unperturbed 
mass dipole: the piece satisfying the same constraint im- 
posed by the gauge condition, dj(M^rii/r 2 ) = 0. A 
good center-of-mass condition might be to set this piece 
to zero at all orders. One can always do so, since one can 
always add a term to that satisfies the same con- 

straint as the unperturbed mass dipole and has arbitrary 
time dependence. 

Further investigation will be necessary to verify that 
this is a good condition to impose. It would not set the 
whole of the corrected mass dipole moment to zero, and 
one should examine the physical content of the nonzero 
portion. To that end, one should examine the physi- 
cal content of the corrections to the monopole moment, 
I?-. , which are symmetric tensors of rank two rather than 
scalars as the mass is, and one should verify how the cor- 
rections to the mass dipole are related to a small transla- 
tion applied I?\/r. Alternative conditions to impose will 
certainly be possible, and at some order, there may not 
be a 'best' choice: any meaningful measure of the 'center 
of energy' will likely cease to coincide with the 'center of 
mass', in which case it would not be possible to set both 
to zero simultaneously. 



D. Effective stress-energy tensor 

Thus far, everything has, for the most part, been con- 
fined to the buffer region. In the region surrounded by 
the buffer, where r ~ m, the equations are not valid 
for multiple reasons: the expansions used are not an ac- 
curate approximation, since the body's field dominates 
over the external background in that region; the region 
B in the external background spacetime may not be dif- 
feomorphic to the region surrounded by the buffer in the 
full spacetime; and unless the body is a black hole, the 
stress-energy tensor of the body will be finite somewhere 
in the region. However, we may ask the question l if the 
solutions are taken to be valid for all r > 0, what stress- 
energy sources them?' This may or may not correspond 
to a small-e expansion of the body's actual stress-energy 
tensor, but it will take a physically meaningful form. 
It will also be relevant in Sec. IIIIEI when defining an 
effective-source equation that may be solved numerically 
for the regular field. 

Now, by construction, the minimal inhomogeneous so- 
lution in the solution to Eq. (|T3|) is sourccd by — 

l^ 7rr (n)' an d the minimal pth solutions for p > are 
homogeneous solutions even at r = 0. This leaves only 
the minimal pth solutions for p < — that is, the parts 
of the singular field that form homogeneous solutions for 
r > 0. For each such I follow the approach taken 

by Gralla and Wald at first order [l2j, defining the dis- 
tributional effective stress-energy tensor 



The right-hand side can be split as 



(62) 



1, \ p n,L / 



(63) 



where the first term is the most singular and 

A (e^" po-(p)4 > (n)) carr i cs au the l ess singular terms. 

First examine the most singular term. Using the 
identities d L r- x = {-l) l {2£ - l)!!^rr and 8%^ = 
—4:TtS 3 (x), we have 



d l d,. 



1 (n) nL _ W l L > 1 (n) 



(2^-1)!! 



d L S 3 (x). (64) 



One can find the covariant form of the right-hand side 
by integrating against a test function ?/> M „. Doing so, we 
have 



^I^ L d L 6 3 (x)V=gd 3 xdt 

= (-I f I d L (V=#^/£f ) dt (65) 
= {-If [ (^C) ai "' ai ) dt, (66) 
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where g = det(g M „), and in going from the second line to 
the third I have utilized the identity di\J—g 



-g- 



Here I have defined I, 



1 to be the tensor that agrees 
with if all en are spatial indices and zero otherwise, 



(n) 



meaning that I,^ x ' 

& (n) 

that /^ ai - a '-««« a , 
shows that 



I^ L d L 6 3 (x) = I I^- a <6(x,z)., ai ... ae dT, (67) 



is STF with respect to g^ v and 
for all 1 < i < I. Equation (l66j) 



where 



^^^(^(r))^^)) (68) 



is the covariant Dirac delta function, and I have made 
use of the fact that t = ron the worldlinc. 
Equation (|63[) now reads 

& (xr(p)9( n )— ^ J 1 ( n ) 0[X,Z). ai ... at OT 



(69) 



Now note that, by construction, A (^-E'' i %o-(p)0(„) J van- 
ishes pointwise for r > 0. If it is nonvanishing as a distri- 
bution, it must have support only on 7, in which case it 
must be proportional to <5 3 (x) or a derivative thereof; but 
from the calculation just performed, that would lead to a 
homogeneous function in the solution. By definition, no 
such functions do appear in the minimal pth solution. So 

A (e^ pa-^^^j vanishes as a distribution. This should 

also be intuitively obvious, since each term in the minimal 
pth solution is constructed directly from its most singular 
term, Ij"^ nt/r^ 1 ; only that most singular term needs 
a source. 

Therefore, Eq. (|()2"j) becomes 



(p) 



(-1)^ 



Since 4^n) con t&ins minimal pth solutions for p = 
— n, 1, the total nth-order stress-energy tensor is 
given by 

T in)=Y,jki^}. 1 i{:r' ae ('r)5(x,z)., ai ... ai dT. 



(71) 

This includes both the moments of the unperturbed body 
and the corrections to them due to 'gravitational energy'. 
If the corrections are added to the moments of the unper- 
turbed body, then one has the total stress-energy tensor 

TliV =£t$ = W I^-<*'6(x,z)., ai ... ae dT, (72) 
n e J t 



where I have defined the normalized corrected moments 

(-iy 
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4(2^-1)!! C«) 



(73) 



It is noteworthy that Eq. ([72]) agrees with the traditional 
form for the multipolar expansion of a material body's 
stress-energy tensor [44T - l46l | . 

Equations (fTTj) and (j?2")) have a physically intuitive 
form: in the limit of small size, a body appears as its 
skeleton of multipole moments (including gravitational 
corrections to them). However, the fields sourced by this 
stress-energy tensor are only meaningful approximations 
in the buffer region and beyond, and it sources only the 
homogeneous (for r > 0) solutions to the wave equation. 
So we can say more precisely that at distances r ^ to, 
the body affects the gravitational field only through its 
multipole structure, and in that sense, it appears as a 
point particle equipped with multipole moments. 



E. Puncture scheme 

In practice, of course, one is interested in obtaining a 
global solution, not merely a local solution in the buffer 
region. Furthermore, solving the field equations in the 
buffer region will not determine the regular field </> Rmi/ ; 
to do so, one must impose boundary conditions, over 
and above the matching condition. To obtain such a 
global solution numerically, and at the same time deter- 
mine the regular field, one can use a puncture scheme. 
Here, I follow the description of Dolan and Barack [30l ]. 
but once a local, analytical expression for the singular 
field is constructed, any puncture scheme may be used to 
numerically find a global solution for the physical field. 



First, define a fcth-ordcr puncture function 



be 



V) 



to 



truncated at order 



(74) 



This function may be obtained analytically by transform- 
ing the buffer-region expression for 



V) 



into convenient 



coordinates. Similarly, define a residual function 



9 [fc](n) 



V) 



9 [fc](n) 



(75) 



With these definitions, ^ryrn) * s a ^* f unc tioii, and it 
agrees with 4\Ja V through order 

Next, define a hollow timclikc tube T of some numer- 
ically reasonable radius around the body. Inside T, one 



9 My notation differs slightly from that of Ref. [30l | , where the label 
[k] corresponded to a puncture containing k orders, which would 



agree with through 



order r 



k — n — l 
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can replace the physical problem, in which one would en- 
counter the material body or black hole, with the effective 
problem 



a P a( P[k]{n) ~ 



167r(T ( My 



(n) -r \n) 



iVpa 



- 16nS [k](r 



(76) 
(77) 



This problem can be solved for <$^° n y By construc- 
tion, the singular terms in E^ v pa^Tk]^,) cancel those in 
Cf \ — l€m(T?K + t, m >). Since the source terms contain 

(n) v. (n) {n)> 

two spatial derivatives of the field, the effective source 
^\k](n) ^ s a C fc " 2 function. So we require k > 2 to ob- 
tain a C° source for the residual field. However, a source 
that is finite but discontinuous will suffice in practice, 
meaning a k = 1 puncture will suffice. 

Outside r, one may solve the physical problem 



(n) 



(78) 



for the physical field </>(^)- When crossing T, one may 
change variables from the residual field to the full field 
via 



V) 



9 [fc](n) 



°[k](n 



(79) 



that is, across T one goes from the effective problem to 
the physical problem by, roughly speaking, 're-inserting' 
the body, adding the analytical expression for 



,Vpa 



to 

the numerical result for (f^Z^ ■ Likewise, when crossing 
in the other direction, into the region bounded by T, one 
may change variables from the full field to the residual 
field by subtracting the puncture. 

Recall that will depend on the lower-order reg- 

ular fields 4\n'<n)' ^° a ^ eacn timestcp in a numerical 
evolution, one must first calculate the first-order residual 
field from the first-order puncture, then use that residual 
field to calculate the second-order effective source, and 
so forth. 

There is obviously great freedom in formulating and 
implementing a puncture scheme. For example, rather 
than taking the truncated expression for <f>^ as a punc- 
ture, one could use any other function that agrees with 
(j),^" through order r k for small r. Additionally, the 
puncture function could be chosen to smoothly go to zero 
at large distances 28 1 or at some finite distance [2£| from 
7. In that case, at sufficient distance the residual func- 
tion would become the physical field, meaning the phys- 
ical field could be obtained by solving Eq. (1751) globally 
rather than only inside V . 



IV. EXPLICIT OUTER SOLUTION IN BUFFER 
REGION THROUGH SECOND ORDER 

I now present the explicit results sufficient to imple- 
ment a second-order puncture scheme. I choose 



pa 
D P 



h^, (80) 
0, (81) 
0. (82) 

The first- and second-order wave equations then read 

0, (83) 

(84) 



F^ v h pa 



The second-order Einstein tensor, <5 2 G A " y , can be written 
with the gauge condition on already imposed or not. 
Here, I use the full expression, (|A3[) . without simplifying 
it via the gauge condition. The solutions to these equa- 
tions are constrained by the first- and second-order gauge 
conditions 

= 0, (85) 
S°^ u h^ =-5 1 V v h%y (86) 

that is, in coordinates centered on 7 the divergence of 
at time t must vanish if a^it) is set to a^(t), and 

the divergence of + e 2 h^ at time t must vanish if 
a p (t) is set to ea^(i) in the first expression and in 
the second. 

In the buffer region, I use Fermi- Walker coordinates, 
in which the spatial coordinates x a span a three-surface 
that intersects 7 orthogonally, and the time coordinate 
is equal to proper time on 7 at the point of intersection. 
More explicitly, at a point x G .Me near 7, the spa- 
tial coordinates are given by x a {x) = — e^,(x')a' a (x, x'), 
where x' is a point on 7 connected to x by a unique 
spatial geodesic (3 that intersects 7 orthogonally, e" is a 
spatial triad orthogonal to u a on 7, and cr(x, x') is one 
half the squared geodesic distance from x' to x. The time 
coordinate is given by t(x) = t(x'), the radial distance 
is given by r(x) = ^/5 a bx a (x)x b (x) = y/2a(x, x'), and 
the radial unit vectors are n a (x) = x a {x)/r(x). Further 
details about the construction of the coordinates can be 
found in Ref. 

Through order r 3 , the metric in Fermi- Walker coordi- 
nates is given by 



9tt = -1 - 2a,ix' ~ (Roioj + a l a j )x l x J 
1 

~ 3 

9ta = --RoiajX'x 3 - - 



(4Roio 3 -o* + iioioj-jfc) arVs* + 0(r 4 ), (87) 
RoiajakX l x 3 x k 



- -Roi aj \ k x z xix k + 0(r 4 ), 
9ab = 5 ab - ^Raibjx'x 3 - ^R albj \ k x l x° x k + 0(r 4 ), 



(89) 
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where the pieces of the Riemann tensor are evaluated on 
the worldline and contracted with members of the tetrad 
(u a ,e%). For example, Ro ia j\k = ^a^;p| 7 " Qe f e f e j e fc- 
An overdot will indicate a covariant derivative along the 
worldline, Roiaj = Ranfiv,p\ u a e^ efe^u p . Because the 
background is Ricci-flat, the components of the Riemann 
tensor and its first derivatives can be written in terms of 
Cartesian STF tensors £ a b, B a b, £abc, and B a bc'- 



£ab — RoaOb, 
B ab 



c pg. 



(aRb)0pqi 



£abc — STF R 0w 
abc 



£abc = g STF e P9 a i?(,0pg|ci 



Ob\c 

pq 



(90) 
(91) 
(92) 

(93) 



where 'STF' denotes the STF combination of the indi- 
cated indices. £ a b and B a b are the even- and odd-parity 
tidal quadrupole moments of the background spacetime 
in the neighbourhood of 7, and £ a bc and B a bc are the 
even- and odd-parity tidal octupolc moments. Identities 
for decomposing each component of the Riemann tensor 
and its derivatives in terms of these tidal moments can 
be found in Appendix D3 of Ref . [47j . 

Following the procedure of the previous sections, in 
terms of these coordinates I assume the trace-reversed 
fields have the expansion 



(n) Z_^i 

p> — n,q,£ 



r p {\nr) q h<t uL 



(t)n L . 



(94) 



An advantage of Fermi- Walker coordinates is that their 
simple geometric properties afford an easy transforma- 
tion to arbitrary coordinate systems. Such transforma- 
tions are available in the literature [48T |. Alternatively, 
one may use the geometric definitions of the Fermi- 
Walker coordinates to express the above expansion in 
terms of the covariant quantity a(x,x'); from there, one 
can express the result in any desired coordinates. There- 
fore, once the singular field is known in Fermi coordi- 
nates, it may be easily transformed into any desired coor- 
dinates to form a puncture function useful for numerical 
implementation. 

In order to obtain a solution for through order 
r, as is necessary for a puncture scheme, I require four 
orders in r in all my expansions, since begins at 

order 1/r 2 . Because the general procedure was shown in 
Sec. Mil and the first three orders in r were derived step 
by step in Rcf. [13] here I shall simply state the results 
for the most part. 



10 Reference [Tj 
and h( 2 ")> 



,0) 



,(2) 



(2) 

and in the case of h^i it did so only after substi- 
tuting the expansion for a M . However, the steps involved in the 
calculation are essentially identical. The results presented here 
were actually obtained by trace-reversing the results reported in 
Ref. [23| ; I present the trace-reversed fields to mesh with the 
foregoing discussion of the nth-order problem. 



First order 



Following the steps outlined in Sec. lIII Al one finds the 
general solution to Eq. (|83j) in the buffer region has the 
form 



(i) ~ "(1,-1,0,0) ^ {\,o,o,i) L ^ (i,i,o,£)"- L 
e=o e=o 



+ r 



(95) 



e=o 



where h^_ lj0 ^ y ^(1,0,0,0)' 



ioi)' and h^-^ 2) are the 



homogeneous modes, each of which defines a minimal pth 
solution to Eq. (|53")) . and of which the other functions 



(i,p,o,.9 



arc linear combinations. 



1. Singular and regular field 



The general solution can be conveniently split into two: 



n (i) 



ft (1) . 



(96) 



Following the definitions introduced in Scc. lIIIBl I define 
the singular field to be the piece of made from a min- 
imal pth solution for p = — 1 — that is, the solution fully 
determined by _ x Q ^ alone. It contains the following 
pieces of the full first-order field (|95l) : 



iSfiu _ 1 rut- 

"(l) " r "(1,-1,0,0) 



"(1,0,0,1)" 8 



(1,1,0,- 



e=o,2 



E^(M,0,^ + O(r 3 ). 



(97) 



(=1,3 



The monopole moment 1^ = h^ l _ 10Q ^ is constrained 
by the gauge condition ([56]) to be of the form I^u^u", 
where 1^ is a constant. The matching condition 

hjlu ^ = 9\^ v then determines 1^ = 4m, where m is 
the ADM mass of the inner background spacetime gi^. 
Given these constraints, Eq. (|97|) explicitly reads 



4m 
r 



h SU _ — _ lQ ma l n . 



28 7 

— did* + -(15aV - 4£ ij )n i: j 
3 6 



1 



19a,£ ij - &V - bQa^^ 



105 



15a £ J — —£ " " " 



a a a \ riijk 



(98) 
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hflf = 2mr (d t a a + B ic e a Jh ij ) 



73 7 

a b viid ai ribd 

30 £ bda B 30 £ dabB 



--e M d a b B ad + —E al - 5a l d t a a - 2a a d t a l ) m 
6 15 



nijk 



hftf = Amr (a a a b - E ab ) + 2mr 

2,, 



£ abi + 2a ( a E b ^ 1 - -B d{a e b)i d 
3 



(99) 

{£ ab - 3a a a b y 

(100) 



n , 



where dtin = 0. 

The regular field is the sum of minimal pth solutions 
for p > appearing in h^s — that is, the part of that 
is determined by global boundary conditions rather than 
the matching condition. It contains the following pieces 
of the full field: 



"(1) ~ "(1,0,0,0) T " t (l,l,0,l)"' J 

+ r 2 ( 

+ V (1,2,0 



2)' H J ^ (1,2,0,0) 



+ 0(r 3 ). (101) 



We can write this as a Taylor series 

> r=0 

1 



"(1) ~ "(i) 



r =o « 



o fyi) <*i 



r=0 



r=0 

xV + <3(r 3 ), 



(102) 



where the coefficients are related to the pieces of the full 
metric as 



7 R/tt/ 
"(1) 


r=0 


- h? v 

(1,0,0,0)' 


(103) 




r=0 


_ IM"> 

~~ "(1,1,0,1)' 


(104) 


TRliv,{ij) 
n (l) 


r=0 


~ z "(l,2,0,2)' 


(105) 


TRfiv i 
"(1) ,i 


r=0 


~~ °"(1, 2,0,0)- 


(106) 



As noted above, all of the functions on the right-hand 
sides must be independently determined by boundary 
conditions, except ^(^2 0)' wrnc h is given in terms of 
the others by 

°"(1,2,0,0) — 2 00 "(l)a& + a a "(l)a& 



1 IRtt 

2 ft (i) 



(107) 



clta _ 3- Rtab R &cn IRd 1 3 a &TRt 

b "(l,2,0,0) - ~2 (!) <^ ft (l)6 +2 W b 

1 



-,&7Ra 1 „aiRtt 1 ^ iRta 
"(1)6, t + a ft (l) ,i + 2 ft (l) > 



gfc 6ft(i) +2"(i) a ,*+2 ft (i) a& '* 



(108) 



-«"(1) ,cO + o f "(l)c+ a a "(1 



R6) 



U ' t (l,2,0,0) 9-"UJ > ti ~ 1 3- 

ccdrabTR cabTRc OR^( fl ^) T"Rtc 

— c n,^ cd -c "(i)c — ^° e cd"(i) 

+ 2a<°^ ) 1 t + ifcf$V (109) 
Here /i^ and its derivatives are evaluated at r = 0. 

3. Stress-energy tensor and puncture scheme 

From Eq. (|7T|) , the body's effective stress-energy tensor 
is determined to be 



= / mu tJ 'u u S(x, z)dr, 
'1 



(110) 



which is the stress-energy of a point mass moving on 7 
in the external background spacetime. Therefore, at first 
order and at distances in the buffer region or greater, the 
body appears as a point mass. 

Following Sec. MI El one may use a puncture scheme 
to obtain the unknown pieces of the field (i.e., h^") in 
the buffer region, along with the global solution ev- 
erywhere outside the buffer region. After transforming it 
into a numerically useful coordinate system and truncat- 
ing it at order r 2 , the singular field given in Eqs. (j9"51) 
POOP defines a puncture hV 2 ?Zy Inside a tube T about 

the body, the residual field ^pf(i) ^ ^(f) ~ ^pf(i) can ^ e 
found by solving the effective-source equation 

puv TlZpa _ la^rp^v _ -puv per (111) 

a P< T ' t [2](l) ~ i07rJ (l) ^ P cr "[2](l) V 111 ) 



°[2](1)' 



(112) 



where ^pui) i s a C° source that can be calculated ana- 
lytically. Outside r, the physical field can be found 
by solving (|83|) . When crossing from one side of T to the 
other, one changes variables between and using 



7 tiv 7 IZpcy 1 tV 'per 

n (l) - n [2](l) +n [2](l)' 



(113) 



Given reasonable initial conditions and a forward evo- 
lution in time, the field obtained with this pro- 
cedure will describe the retarded solution to Eq. (JS3J) 
with the appropriate singularity structure, and on 7 the 
value of the residual field hw^U and its first and sec- 
ond derivatives will agree with those of the regular field 
h^i" ■ Details of the puncture scheme can be found 
in Refs. [H, |30|, HH . Alternative, equivalent first-order 
puncture schemes are described in Refs. (29,. [32l. l33l| . 

Note that because the point-particle stress-energy ten- 
sor (|110[) is a well-defined distribution, the global solution 
for the physical field may be obtained directly from 



E^pji'l^ = -16nT$. 



(114) 
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This field equation was traditionally the starting point 
for approaching the first-order problem. A variety of nu- 
merical methods exist for solving it [Bj], the above punc- 
ture scheme being just one among them. Analytically, 
with arbitrary initial conditions on the Cauchy surface £ 
intersecting 7, the retarded solutioiO reads 

hft =4m f G^^^i/u^'dT + h^ (115) 

where G^ v p i v i is the retarded Green's function for the 
wave operator E^ v pa (with the normalization used in 
Ref. |5|), the integral over 7 is truncated at £, and 
^fi)s ^ s ^ ne contribution from initial data. As was shown 
in Ref. [T^jF 2 ! when one writes this solution in Fermi- 
Walker coordinates and decomposes it into the multi- 
polar form (|95p . consists of tail terms (integrals 
over the past history of the worldlinc and contribu- 
tions from initial data) plus local terms involving the 
acceleration and tidal moments. Furthermore, one finds 
that at least through order r, and modulo any incoming 
waves contained in the initial data, h^ 1 ' is equal to the 

Detweiler- Whiting regular field 4m J G Ktlv ^^u^ u u dr. 
where G RMI/ M ' y ' is the Detweiler- Whiting regular Green's 
function (2lj . Correspondingly, through order r, h^i" is 
equal to the Detweiler- Whiting singular field. 



3. Equation of motion 

Evaluating Eq. (|6Tj). one finds that the zeroth-order 
term in the acceleration is 

< ) = °- ( 116 ) 

Therefore, a M = 0(e), and at leading order the body 
behaves as a test mass, satisfying the geodesic equation 
in the external background spacetime. Corrections to the 
geodesic equation appear only at subleading order. 



B. Second order 

Again following the steps outlined in Sec. IIII A[ one 
finds the general solution to Eq. ([84]) in the buffer region 



11 By 'retarded solution', I mean the causal solution obtained with 
fields propagating forward in time from arbitrary initial data; I 
do not mean the particular causal solution corresponding to zero 
incoming radiation at the infinite past. Given that the influ- 
ence of initial data decays with time, any causal solution should 
approach the latter. 

12 In Ref. [l3ll . the equality was shown only for a M = 0; it is shown 
for arbitrary acceleration in Ref. [EJ. 



has the form 

1 2 1 3 _ 

"(2) ~ r 2 Zv (2,-2,0/)" r Zv (2,-1,0/)'' 

1=0 1=0 

4 5 

+ X] ^(2,0,0,£)"- L + T Z\2t ^(2,1,0/) 
1=0 1=0 

+ hr[%>,i,o) + r %i,i) n <] + °( r2 )' 

wneie r^ 2 ,-2,0,i) ~ i (2) > "(2,-1,0,0) — J (2)> "(2,0,0,0)' anu 
^(21 1) are *^ ie homogeneous modes. The other func- 
tions e) are hnear combinations of those homoge- 
neous modes plus quadratic combinations of the homo- 
geneous modes appearing in . 

1. Singular and regular field 

Following Sec. IIIIB1 I split the second-order general 
solution as 

h% + (118) 

The singular field is the sum of three terms: 

*\2) = (-2)^2) + (-1)^(2) + h™)", (119) 

where each of the terms is given through order r in Ap- 
pendix [D] (-2)^(2) an< ^ (-1) (2) are ^ e mm i ma l pth so- 
lutions for p = — 2 and — 1, and h(p\ V is the minimal 
inhomogencous solution to 

E^pahf^ = -2 (s 2 G^[h {1) ] - PG^frfv]) . (120) 

In other words, h,^ is the sum of all terms in that 
directly involve the moments I(2)i ano - ^(2V ^^ 1C 
regular field h,£" is the sum of all minimal pth solutions 
for p > 0, plus the minimal inhomogeneous solution to 

E^ P ahf 2 p } a = -2S 2 G^ [hfc]. (121) 

That is, h^2)' contains all terms in that arc not 
determined without imposing boundary conditions be- 
yond the matching condition. The regular field hf 1 ^ 1 ' = 
eh^ v + e 2 h^ v is a free gravitational field, a smooth so- 
lution to the homogeneous EFE through second order in 
e. 

With these definitions, the singular field contains 
the following pieces of the general solution: 
2 3 

TSpv _ J_ V bT L fl L + - V h >tVL h L 

"(2) ~ r 2 (2-2,0, i) 11 + » Zy (2,-1,0,/)" 

e=o 1=0 

4 

Z^ (2,0,0/)" ^ 2—1 (2,1,0/)" 
t=l £=0,2,3,4,5 
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First examine ^_ 2 ) ^(2) ' which is the homogeneous so- 
lution of the form I^n.i/r 2 + 0(r _1 ), defined to con- 
tain only those terms in NL\ that are directly pro- 
portional to the dipole moment I^) ■ After splitting 
the dipole moment into its mass and current pieces, 
1(2) = ^(2) + ^(2) 1 ' tnc S au S c condition (|56| yields 
M^h^ = 0, which constrains the mass dipole to have 



the form 



'(2) — •* J "(2) "t- 3 u " u 4j J "(2)- g au g e 
refinement can be used to eliminate the ta component 
of and the ab component of [l]| , leaving the 

dipole moment in the canonical form 

Au' l u v M i + Au^e v ^e aji Sj, (123) 



(2) 



where I have defined M l = M%l/4 and Sj = e 3 -«,S%/4. 
The factors of 4 are chosen based on the matching condi- 

(2 —2) (2) 

tion h^v = g\', which determines S l to be the ADM 
angular momentum, and M l the standard mass dipole, 
of the inner background spacetime gi^v I set M l = to 
enforce the condition that the body is effectively centered 
on 7, which leads to 



(-2)^(2) 



5 hO(r i ). 



(124) 



The gauge condition (|56|) further determines that the 
spin is constant; dtS 1 = 0. 

Next, the monopole field (-i)h^ is the homogeneous 
solution containing only those terms in that are di- 
rectly proportional to the monopole moment 1]%., which 
I now denote 8m tiv . It has the form 

■1)^ = ^+0(1), 



(-l)' l (2) 

where the gauge condition constrains Sm^" to be 
1 



dm 



g 7Tm (l) 1 
2 T T?«^ 2 



(125) 

(126) 
(127) 



-m8 ab hfjf. (128) 
3 



Sm ab = -mfcfff - -m5 ab h^ 
3 (1) 3 (1) 

This is a gravitationally induced monopole; it can be 
thought of as the perturbation of the smooth field hJt" 

by the body's local field L*!", or as a correction to the 
body's local field h,^ by the smooth field hJt u . Mathe- 
matically, it arises because fya/*" requires the addition of 
a 1/r monopole in order to satisfy the gauge condition. 
One could add an arbitrary, time-independent monopole 
of the form u^u v 8m — that is, one that behaves exactly 
like the first-order mass — but I choose to incorporate that 
term into the definition of m. 

The regular field contains the remaining pieces of the 
general solution: 

KzT = ^0,0,0) + ^S,o,i)^ + 0(r 2 ); (129) 



terms quadratic in h^t" , arising from the source in 
Eq. (|12ip . would appear at order r 2 . The functions in 



Eq. (|129p are constrained by the gauge condition 
which relate them to each other and to the first-order 
regular field, but they are otherwise undetermined un- 
til further boundary conditions are imposed. As at first 
order, we can write this as the Taylor scries 



, R/J.V _ 1 Kfiv 
"(2) — "(2) 



r=0 



, R/J.V 
n (l) , 



r=0 



x l + 0(r 2 ), (130) 



where the coefficients are related to the pieces of the full 
metric perturbation as 



rltftu 
n (2) 

n (2) 



r=Q ""(2.0,0,0)' 
„ ~ "(2,1,0,1) " 



(131) 
(132) 



2. Stress- energy tensor and puncture scheme 

From Eq. (fTTj) the body's second-order effective stress- 
energy tensor is 



1 (2) 



u^S v)a 5{x,z), a dT + / ^-5m ,lv 5{x,z)dT, 



(133) 



where I have defined S a ^ = e"e^e abl Si (not to be con- 
fused with ). The first term is the effective stress- 
energy due to the body's spin; the second, that due to 
the induced monopole Sm^ . If the spin vanishes, then 
at distances in the buffer region or greater, the body still 
appears as a monopolar point mass, but with a nonlinear 
gravitational correction to the monopole field. 

Again, a puncture scheme can be used to obtain the 
unknown pieces of the field in the buffer region, h^ ) 

and \ along with the physical field globally. The 
singular field given in Appendix (|D|) . truncated at order 
r, defines a puncture hm^y Inside the tube T, the resid- 



ual field h 



TZpa 



"Ppa 



*[i](2) = "(2) — "[i](2) can ^ e found by solving 
the effective-source equation 



a P°' l [l}(2) 



-l&-KT^ ) - 25 2 G^[h (1 y) 

-^*V% (134) 

= s [i](2y ( 135 ) 

where is a discontinuous but finite source that can 

be calculated analytically. The analytical expression for 
the source will involve the first-order residual field, which 
is to be determined numerically at each time step along 
with the second-order residual field. Outside T, the phys- 
ical field can be found by solving (|54"|) . When cross- 
ing r, one may change variables between h^) and 
using 



n (2) 



illpa 
n [l](2) 



,Vpa 
n [l](2)- 



(136) 
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This scheme is essentially identical to that at first or- 
der. In principle, it can be implemented almost imme- 
diately in existing codes. One need only transform the 
puncture function into a convenient coordinate system 
and then calculate an analytical expression for the effec- 
tive source in those coordinates. 

Unlike at first order, where one could write down an 
analytical solution in terms of a retarded Green's func- 
tion, at second order no analytical solution presents it- 
self. The part of the retarded solution that is sourced 
by 5 2 G p " l/ cannot be written as an integral over all space, 
of the form J G pv P > IJ >5 2 G P a dV, since the 1/r 4 singular- 
ity in <5 2 G A " y would make that integral diverge at every 
point in space. The actual solution consists of a volume 
integral over f2 plus a surface integral over dfl fl3j ] : 

an 

" 27 X Ga0 t' s ' 52G1 ' s ' M dV '> ( 137 ) 

which requires one to know the field on the timelike part 
of dfl in the buffer region. Therefore this formal solution 
is not especially useful. However, one could use to it to 
find analytical expressions for the homogeneous modes 
that appear in h,£ v . Using the analytically determined 
solution in the buffer region as data on the part of dfl 
embedded therein, the integrals could be expanded for 
small r with methods from Ref. [lj]; decomposing the 
result into STF pieces, one could then read off the homo- 
geneous modes. This small-r expansion is complicated 
by the fact that an expansion of the Green's function 
would rely on the Hadamard decomposition [f|, which is 
valid only in a convex normal neighbourhood, while the 
volume integral in Eq. (|137[) gcnerically extends over a 
much larger region. So, to evaluate the volume integral 
at a point x, one must restrict one's attention to points 
x very near the initial time-slice E, such that the volume 
in which G a ^ 7 'gi has support inside the integral lies en- 
tirely within the convex normal neighbourhood of x. I 
leave this calculation to future work. 

However, the part of the retarded field that is sourced 
by can be immediately written down. It reads 



agree with 4 f S7 a >(G R ^ p , a ,uP' S a ' a ')dT at least through 
order r, but I have not performed the explicit expansion. 

3. Equation of motion 

Evaluating Eq. (pT|) and then making use of the gauge 
condition 8 G {d t h^ + 9f>/im) | r _ = 0, one finds the first- 
order term in the acceleration is 

ma* (1) = ^m5°(5abhffi,o,i) + + mSPdthfa - B ij S 5 . 

(139) 

In terms of the regular field hf^ , the result reads 



'(i) 



(140) 

This is the usual first-order result. The first three terms 
are the first-order gravitational self-force, and the last 
term is the Mathisson-Papapetrou spin force. 



Summing the terms in the acceleration as a p 



ea^i + 0(e 2 ) and converting from h R)J - 1 ' to h^", one ar- 



"(i) 
rives at 



;(g^ + u»un(2h*. ( ,-h%. v )u<>u° 



+ —R^p^S^ + 0(e 2 ). 
2m 



(141) 



Through order e and excluding the spin, Eq. (|14ip is the 
geodesic equation in the effective metric g^ v + h R v @|. 
Through that order and with the spin, it is the equation 
of motion of a test body in g^ v + h^ v . In the punc- 
ture scheme described in the preceding section, it should 
be used (after dropping the "Q(e 2 )" ) to sclf-consistently 
evolve the position of the bodvPl 

In going from Eq. (|140[) to Eq. (|141[) . I have replaced 

eh^y with h R u and dropped the 5°. Doing so does not 
alter the acceleration at order e. However, recall that 6° 
indicates evaluation at a M = 0. This evaluation is to be 
performed only at time t, leaving the past history of 7 
unchanged. If ft^" is evaluated analytically, in terms of 



4V a <(G^ 



p' a'l 



S°<*) + G^ plc ,8m pa dr. (138) 



This can be expanded for small r, decomposed into mul- 
tipoles, and its contribution to hf"^ determined. (Its 

contribution to is simply the minimal pth solu- 

tions for p = — 2 and p = — 1.) At least through or- 
der r, the 8m} LV term in h,£ is found to agree with 

/ G^" p > (7 >5m p '' T ' dr. Explicitly, this consists of (i) tails 
of tails, since 8m pc7 is made up of tails, and (ii) tails 
multiplied by local factors of the tidal moments and ac- 
celeration. The spin term in can be expected to 



13 In fact, to be consistent with my equations for the second-order 
singular field, one must use an equation of motion at least as 
accurate as Eq. 11411 1 . If instead the less accurate equation = 
is used, then the mass dipole terms in the singular field cannot 
be set to zero, since if a worldline 7 of vanishing acceleration is 
chosen, then the mass dipole evolves according to 1 1 3H 

^-M" = B?- vp vU v u<>M° + —R p - vpa u v S'> a 



(142) 



If one is interested in finding a solution to the second-order EFE 
on a very short timescale, one can use = 0, include the 
M 1 terms in the singular field (presented through order r° in 
Ref. [H), and evolve M i via Eq. l fl42l . 
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tails, then -\{g^ v + u»u")(2h* - h^ tJ . u )u p u' y includes 
tail terms, which depend on the past history, plus local 
terms involving the acceleration. In particular, there is 
a term — -^-ma M which is the traditional gravitational 
anti-damping term found by Havas (49| (as corrected in 
Ref [13] )■ It might seem that going from Eq. (|140|) to 
Eq. (|141|) has thus transformed a second-order ODE into 
a third-order, non-physical one. However, in practice 
no third-order ODE will arise; the regular field, not the 
tail field, will be found from a puncture scheme (or the 
various other numerical regularization schemes used at 
first order), and Eq. (|141[) will be evolved directly, as a 
second-order ODE. Evidence from the scalar field case 
suggests that no problems will arise when performing 
a self-consistent evolution in this manner: when writ- 
ten in terms of tails, the scalar self-force involves 
terms that would lead to runaway evolutions Q , but self- 
consistent numerical evolutions, using an equation of the 
form (|14ip . show no evidence of runaway behavior [5lj . 
Alternatively, it should be possible to numerically evolve 
= ea^ directly, with given by Eq. P^|) or (jTlUl) . 
Finally, it is important to note that the acceleration in 
the form (| 141[) would have followed directly from the 
gauge condition (|2T)|) . Therefore, assuming Eq. (|141[) 
is well behaved, in the sense that the at 1 term that is 
hidden within the right-hand side gives rise to no non- 
physical behavior, then the assumed expansion (|16[) of 
the acceleration on 7 is unnecessary: at least through 
the orders explicitly examined here, the assumptions on 
the field alone suffice to arrive at a well-behaved equation 
of motion. These issues will be discussed further in the 
sequel (24|. 



V. DISCUSSION 

The analytical calculation of the second-order gravita- 
tional field outside a small, compact body is now com- 
plete, in the sense that sufficient information has been 
found to numerically obtain a global solution outside the 
body. I have defined a split of the second-order field 
into singular and regular pieces, and I have found an ex- 
plicit expression for the singular field through order r 
in a local expansion. In an earlier work [l3j . to derive 
the first-order equation of motion I found the second- 
order singular field through order r°, but that is insuffi- 
cient to find a global solution, and I also dropped terms 
with explicit acceleration dependence. The full results 



This is easily found by using the explicit results for h^ a (r = 0) 
and di/h^lr = 0) in Fermi coordinates in Ref. |14| (also repro- 
duced in Ref. Q). But note that Table I in RefHll and Table 
II in Ref. 0] are missing a factor of 4 from the ma a term in the 
quantity ci 1,0) = hf a {r = 0). The missing 4 appears, correctly 
situated, in Eq. (E.9) from the former reference and (23.10) from 
the latter. 



through order r, first reported in Ref. [23| and made ex- 
plicit here, do suffice to obtain a global solution via a 
puncture scheme. Reference [23[ also showed that the 
regular field, as I have defined it, is entirely responsible 
for the second-order force; further details of that deriva- 
tion will be presented in the sequel [HJ. 

Note that I could have bypassed the split into singular 
and regular fields entirely, simply defining a puncture and 
a residual field. From the perspective of the equation of 
motion, the only matter of practical importance is that 
through order r, the residual field contains every term in 
the field that is smooth at r — 0; that is, given the field 
in the form (|23[) , one can simply take the nth-order punc- 
ture to be everything but fi^o 0) ~^ r( ^(2i i) ni + ^( r )> 
where the choice of the 0(r 2 ) terms has no practical im- 
pact. However, the definition of the singular and regular 
fields through all orders in r affords a prettier picture, in 
that the regular field, as I have defined it, is a smooth 
solution to the vacuum Einstein equation, and the effec- 
tive metric g^ + h^ v in which the motion is geodesic is 
therefore a physically meaningful spacetime. 

More generally, I have shown how the same procedure 
works at arbitrary order in a large class of gauges. I 
defined an nth-order split of the field around a small 
body into a singular/self- field and a regular/effective 
field, where the singular field involves the body's mul- 
tipole moments and corrections thereto, and the regular 
field is a smooth solution to the vacuum Einstein equa- 
tion, which is determined by global boundary conditions. 
I also described how to use the former to find the latter 
(and the global solution) with a puncture scheme. In one 
sense, this split into self-field and effective field marks an 
improvement over the non-pcrturbativc split defined by 
Harte [lH, since my regular field satisfies the vacuum 
EFE at all orders in e, while his does not. However, it 
lacks the elegance and covariant utility of his definition, 
since it is defined in terms of a series solution in a par- 
ticular coordinate system. 

I also found the form of the body's nth-order stress- 
energy tensor as it appears (in a sense) from distances 
in the buffer region and beyond. And I have shown how 
an nth-order equation of motion follows from the nth- 
order gauge condition, though an appropriate constraint 
on the 'corrected dipole moment' must be imposed to 
make this a meaningful and useful representation of the 
body's motion. 

From these considerations, one sees that a specifica- 
tion of the body's multipolc moments (including the 
aforementioned dipole moments and together with ini- 
tial data) is sufficient to find the physical metric every- 
where outside the body and to self-consistently evolve 
the body's position. 

A. Comparison with other work 

Besides Refs. [TH, [23| and the present paper, the main 
work done on the second-order self-force problem has 
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been performed by Rosenthal [§3, HH , Detweiler (3(| , and 
Gralla [23 |, The most immediate, though relatively mi- 
nor, difference between those works and the present pa- 
per is that they assume the small body is a monopole 
through second order, while the present work is valid 
for a spinning body. However, all of these studies share 
common aspects: gathering information about the sin- 
gular behavior of the metric perturbation by matching 
to an inner expansion; and expressing the force in terms 
of some specified, regular part of the field. Here I will 
discuss only the methods of finding the second-order per- 
turbation, leaving discussion of the force itself to the se- 
quel 0. 

Rosenthal found an analytic, global (excluding 7 it- 
self) solution to the second-order Einstein equation [3J], 
along with the general form of its divergent pieces near 
7 [35| . Since he focused on finding an expression for the 
self-force, he did not derive an explicit local expansion 
for a singular field that could be used within a numeri- 
cal puncture scheme. However, the larger disadvantage 
of his approach is that his solution, and his method of 
solution, relics on a peculiar choice of gauge in which 
hn\ V ancl dpt\i) V both vanish on 7. This means, most 
strikingly, that the first-order self-force vanishes. (It also 
means that Sm til ' vanishes.) In such a gauge, there is no 
obvious way to calculate the metric perturbation in prac- 
tice. Also, the perturbation will grow large with time. 

Recently, Detweiler has taken the opposite focus, not 
deriving an expression for the force but instead discussing 
general features of a second-order split into singular and 
regular fields [3(|. Like the present paper, he advocates 
using a locally obtained singular field together with a 
puncture scheme^ His discussion of the second-order 
problem has the advantage of being gauge-independent, 
though it docs not provide an expression for a second- 
order singular field or an algorithm for constructing one. 
If his singular field is to be taken from the metric of a 
tidally perturbed black hole, as indicated in Appendix 
A of Ref. J36|, then the metric perturbation will be in 
an inconvenient gauge. When transforming to another 
gauge, one must then decide how to split the transformed 
perturbation into singular and regular pieces. Depending 
on how that is done, the motion may or may not be 
geodesic in g^ iv + h^ v ; for example, with Gralla's choice 
of h^ v and h^ v , the motion is not quite geodesic in g^ v + 
[1^. Of course, any choice of singular and regular 
fields will do, so long as one has a convenient equation of 
motion in terms of that particular choice. 

Besides these differences in scope, one detail of De- 
tweiler's approach appears to be at odds with the re- 



5 The equation he presents for the second-order retarded field is 
purely formal, since it includes 5 2 G A " y [h^] at all points down 
to r = 0; because S 2 G ,1V [h^] ~ this has no clear distri- 

butional meaning. But the essential aspect is his effective-source 
equation for the regular field, which is well defined and takes the 
same general form as that presented here. 



suits presented here. Based on the fact that the particle 
obeys the geodesic equation in g^ u + eh^u , he assumes 
that through second order, the stress-energy tensor of the 
body is that of a point particle in g^ + eh^i, : 

T^[h^} = J m«V S A (x a -z a ) ^ (M3) 
7 Jdet [g pcr + e/ipi 1 ^ 

This expression for the stress-energy does not agree with 
my result for eT?X + e 2 T^. One can see that the results 
differ simply by noting that Eq. (11431) contains no com- 
ponents orthogonal to the worldlineP^I while Eqs. (|133[) . 
(|127|) . and ()128|) show that does contain such com- 
ponents. The source of the difference is not clear. It 
might be a quirk in my choice of gauge; the components 
of the regular field, which determine the form of 8m} 11 ' 
and therefore of TftX, can be adjusted via a smooth gauge 
transformation. However, it is not obvious that Eq. (|143j) 
should be the correct choice. If a metric perturbation 
h^y is taken to naively satisfy the ill-defined second-order 
EFE 8G» v [h] = 8nT^[h] - S 2 G tny [h], then geodesic mo- 
tion in g^ v + /i M „ at linear order in h^ u formally follows 
from the Bianchi identity = V 1/ (8ttT^ 1/ [/i] - S 2 G^[h]), 
as shown in, e.g, Sec. IIB of Ref. [l^|; but if the met- 
ric perturbation is taken to naively satisfy the (also ill- 
defined) EFE 8G^[h] = %TrT^[h K ] - 8 2 G^[h], then 
geodesic motion in g^ + h^ v does not so follow from 
the Bianchi identity = V u (8nT^[h R ] - 8 2 G^[h\). Of 
course, neither of these Bianchi identities are well defined 
on 7. So I leave the issue open. 

Most recently, within the context of a power series ex- 
pansion rather than a self-consistent one, Gralla [22[ has 
obtained results very similar to my own: an expression 
for a singular field accurate through order r in a local co- 
ordinate system (centered on a geodesic 70, rather than 
the 7 of the present paper) a prescription for obtaining 
the regular remainder of the field via a puncture scheme, 
and a second-order equation of motion written in terms of 
that regular field. His method of finding the second-order 
field (as well as the equation of motion) begins by finding 
an inner expansion in a gauge that is mass-centered on 
the geodesic 70. That is sensible within the context of the 
short timescales of Gralla's approach, because on those 
timescales, the body's deviation from 70 is of order e or 
smaller, meaning the body can be translated to sit atop 



Note that Detweiler advocates using an ordinary power series 
expansion rather than a self-consistent one, meaning Eq. 111431 

is to be expanded not only for small h^}\ but also for — > 
2 (0) ez (l)' wnere z (o) ls ^ e geodesic initially tangential to z' 1 . 
To compare my stress-energy tensor to that expanded form, one 
would similarly expand tT^ + e 2 T^ y Section IIB of Ref. [3 
contains details of such an expansion, the resulting second-order 
EFE, and its formal solution, in the case where /ij^i is replaced 
by h ( p} in Eq. dT43l ■ 
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70 via a gauge transformation. The solution in that mass- 
centered gauge, though it applies to any spherical body, is 
equivalent to typical expressions for a tidally perturbed 
Schwarzschild black hole, such as those in Ref. it 
is also similar to the gauge used by Rosenthal. After 
converting this inner expansion into an outer expansion, 
the solution in any arbitrary, smoothly related gauge is 
written as the sum of the original solution plus terms in- 
volving the arbitrary gauge vector. That gauge vector is 
implicitly determined, up to a free choice of some initial 
values, by one's choice of gauge for the regular field, and 
it is found by solving a set of simple transport equations 
on 7o. As such, Gralla's second-order field is more gen- 
eral than that presented here, since it is expressed in a 
class of gauges rather than a single gauge. 

The more obvious difference is that his approximation 
is inherently limited to short times. To determine long- 
term effects it will need to be combined with a scheme 
such as the two-timescale expansion of Hindcrcr and 
Flanagan In Gralla's second-order field this limi- 

tation manifests as secularly growing mass dipole terms 
representing the body's movement away from 70. These 
are the terms involving A 1 (corresponding to M z in my 
notation) in his Eqs. (B4) and (B6); in my approach, 7 is 
chosen to ensure that these terms vanish. Less notably, 
as mentioned above, his choice of singular and regular 
fields differs from my own and leads to an equation of 
motion that is not quite the expanded geodesic equation 
in + ■ To determine whether his results agree with 
mine, one would have to find a gauge transformation from 
the Lorenz gauge to his class of gauges and perform an 
expansion of the worldline 7 about a geodesic 70. Alter- 
natively, one could avoid the expansion of the worldline 
by setting a M to zero in all my expressions and reinsert- 
ing the mass dipole terms that were set to zero (through 
order r°, these terms can be found in Ref. [HI). Doing 
so would immediately yield the metric centered on 70 . In 
either case, given the very similar underlying methods, 
disagreement is unlikely. 



B. Future work 

At second order, in the case of a compact body with 
slow internal dynamics, the present paper (together with 
its sequel [24]]) and that of Gralla essentially represent a 
complete solution to the analytical portion of the prob- 
lem. The only obvious analytical work remaining in my 
treatment is to find a closed form expression for the sin- 
gular and regular fields (or for two fields that agree with 
them through order r) — in other words, to find second- 
order equivalents to the Detweiler- Whiting fields. One 
might do this by finding suitable quadratic combinations 
of first-order fields, as in Ref. |34j, to make up a singular 
field. As mentioned in Sec. IIVB 21 to find a closed form 
for the regular field, one could write the retarded solu- 
tion in integral form outside a worldtube embedded in 
the buffer region and then devise a decomposition of the 



original integrals into two pieces, such that one of them, 
when expanded for small r, contributes to the O(r ) and 
0{r) pieces of and the other does not. Details of 

such integral representations and their expansion are con- 
tained in Ref. 

Closed-form expressions might allow numerical alter- 
natives to the puncture scheme. And they would aid 
in comparing my results with other analytical methods, 
such as those of Harte and Galley JEl, once those 
other methods are applied at second order. However, 
this is not strictly necessary. So at second order there 
is only one major remaining goal: to numerically imple- 
ment a puncture scheme. Using the results of this paper 
alone, one can implement such a scheme by simultane- 
ously solving Eqs. (jllip and (|134[) for the regular field 
in a region effectively covering the body (in the sense 
described in the introduction), Eqs. ([53")) and (jM)) for 
the full field everywhere else, and Eq. ()141[) for the po- 
sition of the body. Doing so would yield a solution to 
the EFE through second order, since the relaxed EFE 
([8]) and the gauge condition ([7]) would be satisfied up to 
third-order errors. However, to obtain a solution accu- 
rate on longer timescales, a preferable prescription would 
be to replace (|141|) with a second-order-accurate equa- 
tion of motion. For the case of a non-spinning body, that 
equation is given by Eq. (17) in Ref. (23[. Further details 
of its derivation will be presented in the sequel [24| . 

At nth order, much work remains. Alternative gauges 
could be explored. A different choice of field variable 
(j)^ and gauge terms Z^ up and might simplify the 
equations. In particular, one might find a combination 
of Z M „p, and to eliminate logarithmic terms from 
the solution, as in post- Newtonian theory [4(|. Among 
other benefits, this would allow one to construct a strict 
power series expansion in a generalized wave gauge. 

One should be able to explicitly relate the multipolc 
moments I^ vL (including the corrections to the body's 
unperturbed moments) in the outer expansion to prop- 
erties of the matter and gravity inside the body, analo- 
gous to how radiative multipole moments are related to 
the source moments in post- Newtonian theory [37j ■ This 
might be most easily done in the case of a material body, 
where the moments could be expressed as integrals over 
the body's interior. In particular, one would like to en- 
sure that the corrected mass dipole provides a meaningful 
measure of the center of energy relative to the worldline 
7. It may be helpful to relate the corrected mass dipole 
to a derivative of some linear momentum, defined as an 



integral over a surface around the body [54 1. 

Mass dipole aside, concrete relationships between the 
multipolc moments and properties of the body would ren- 
der their meaning more transparent and aid in formulat- 
ing physically realistic models for their magnitude and 
evolution. Models for both the spin and quadrupole mo- 
ments will be of use in modeling binaries, since those mo- 
ments appear in the equations of motion at the same or- 
der as the first- and second-order self-force. See Ref. [HH 
for a recent discussion of modeling the quadrupole mo- 
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ment of the small body in an EMRI. 

One could also consider a wider class of small bod- 
ies. The assumed inner expansion, in which all lengths 
are scaled by e and time is unsealed relative to the exter- 
nal scales, restricts the approximation scheme to compact 
bodies with slow internal dynamics. The single scaling of 
lengths means that the body's linear size d is of the same 
order as its mass. This makes its multipolc moments 
scale as m e+1 , as opposed to the more natural scaling 
md e . From the perspective of the outer expansion, one 
could simply let the multipole moments scale with in- 
dependent powers of e to describe a more generic body. 
However, there would still be an underlying limitation in 
the approximation scheme, since it assumes the buffer re- 
gion is outside the body; if the body is very diffuse, such 
that d > TZ, then the scheme cannot be expected to be 
accurate. Not scaling time in the inner expansion means 
that no functions in it depend on the fast timescale t/e 
in the way they depend on the short lengthscale r/e. Al- 
lowing fast internal dynamics in the body would result 
in rapid oscillations in the metric. Removing this restric- 
tion would likely require a more intensive and delicate 
approach to the outer expansion. 

Finally, the outer expansion cannot be expected to be 
valid on an infinite domain (though it may turn out to 
be). Since in the end we seek the waveforms produced 
by the body, we may need to match the outer expansion 
to an outgoing-wave solution at infinity. 
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Appendix A: Second-order curvature tensors 

The second-order Ricci tensor S 2 R a p corresponding to 
a metric perturbation h^is given by 

S 2 R a /3 = —^h^-v (2/i M ( Q;/ 3) — h a p.^) + \h^ v -ah^.p 

— \h^ v {2h^ a .p) u — hoip-^v — h^ v . a p) . (Al) 

The second-order Einstein tensor with indices up can be 
obtained using 



(A2) 



+ 9^9 



(A3) 



where I have specialized to a Ricci-flat background and 
utilized the identity 5{g^ v ) = -h^ v . S 2 G alB can be writ- 
ten as a functional of two arguments, 6 2 G a "[h, h], such 



that, for example, 
S 2 G a/3 [h s + h R ,h s 



S 2 G a ?[h s 1 h s ] + S 2 G a P[h s ,h R ] 
+ S 2 G ap [h R ,h s } + S 2 G a ^[h R ,h R } 



(A4) 



With an abuse of notation, I write S 2 G a ^[h] = 
S 2 G a "[h,h] when both arguments are the same. With 
a greater abuse of notation, I write 5 2 G a ^[4)] to de- 
note 8 2 G a "[h [(/)]], meaning h^ u is replaced by its expres- 
sion in terms of M „. For example, S 2 G al3 [h] denotes 
the right-hand side of Eq. (|A3[) with h^ v replaced by 

Analogously, S n G a ^ denotes the piece of G"^^,, + 
hftv] that contains n factors of and its derivatives, 
and I write it as S n G al5 [h, h] or S n G afi [h}. 



Appendix B: STF decompositions 

I here reproduce standard formulas from Refs. |39ll4ll|: 
Eq. (|B18|) is the only identity not contained therein. 

Any Cartesian tensor T s (8, <fi) on a sphere can be ex- 
panded along STF combinations of unit vectors as 



T 



S/ 



oo 



S(L) 



til, 



where the coefficients are given by 

T S(L) = W+W f T S (6 ^ 

Anil 



(Bl) 



(B2) 



where x\\ = x(x — 2) ■ ■ • 1. 

These coefficients can then be put into an irreducible 
form. For example, for s = 1, we have 



rpa{L) _ rfiaL , ja(i e rpL—1) sa(i t rf^L-l) 

where the 's are STF tensors given by 



T, 



L+l 

(+) 



T, 



L _ 
(0) 



rppq(L-l i t ) 



rfiL-l _ 
± {-) ~ 



l+l 

2*-1 t ,-l-i 
21+1 J 



(B3) 

(B4) 
(B5) 

(B6) 



Similarly, for a symmetric tensor T$ with s = 2, we have 



T. 



ab(L) 



S ab K L +f^ 

+ STF STF ( e p a i .r b ( + 1 l 1 + Sai t f$_i 



+ <5 aJf e p ^_ 1 Tp L l ) 2 + 8 ail 8 bil _ 1 f K L J' 2 , \, (B7) 



(-2) 



where 
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(L+2), 



(B8) 
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T, 



(+i) 



21 



L+l 



(0) 



£ + 2 L+l 

6£(2£- 1) 
(f + l)(2£ + 3) ~2 



(B9) 



STFCfy^L-x), (BIO) 



T ^ = (ZTTK27TT) S l T ? (t <^ i - 2e ^- 1 } 



f-(-2) _ 2 ^~ 3 T jib 
1 L-2 = 2i+l i <jfe) L~2 



(Bll) 
(B12) 
(B13) 



These decompositions are equivalent to the formulas for 
addition of angular momenta, J = S + L, which results 
in terms with angular momentum I — s < j < £ + s; the 
superscript labels (±n) in these formulas indicate by how 
much each term's angular momentum differs from £. 

When manipulating quantities of the form (jBip . the 
following identities are useful: 



n c n L = n cL + ^—■ 0^11^-1} , 



(B14) 
(B15) 



From the latter, one finds that til is an eigenvector of 
the flat-space Laplacian, satisfying 

r 2 d l d t n L = -£(£ + l)n L . (B16) 

Equation (|B2j) is related to the integral identity 

(B17) 



f n s n L dn = STF S{nl2 '" 6l ° +e - ll ° +e} 



4~ 



(« + * + !)!! 



if s + £ is even and zero otherwise. Here the curly braces 
indicate the smallest symmetric (unnormalized) combi- 
nation of the enclosed indices. With a bit of combina- 
torics one can use this identity to show 
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'47T 



s\A il i \.. ip i P L B 1 



A S B L i- j n s n L dn = '""^J- (B18) 



(s - £)\(s + £ + 1)\V 



for a symmetric tensor As and an STF tensor Bl sat- 
isfying £ = s — 2p, where p is a nonnegative integer. If 
£ 7^ s — 2p, then the left-hand side vanishes. 



Appendix C: Normalization of logarithms 

In my ansatz (|23[) for the field, I write the logarithms 
with no normalization in their argument. In actuality, 
they must read In — , where r$ must be determined by 

I 



boundary conditions. At least at low orders, the nor- 
malization can be expected to be r = 2em, which is 
what would arise from the effect of the mass on the 
shape of the light cones (i.e., the shift from u = t — r 
to u = t — r* as a lightcone coordinate). However, 
within the present context, the actual choice is arbi- 
trary. Consider the case when Cv\ contains no log- 
arithms. Then a logarithm first arises in the solution 
when the source for the Poisson equation (|27|) contains 
a mode with £ = £ p ; all other logarithms will either 
arise in the same way, or be sourced by logarithms that 
arose in that way. So we may consider those logarithms 
alone. Suppose we begin with a normalization r but 
then change to r' . Then r p h[ Jp \n('i"/ro)4 > n'pi£ * s replaced 
with rPn Lp Mr/rtttfSZ + r p n Lp Hr'Jr,)^. The 
term hi^g/ro^^f/ 3 is then simply absorbed into the ho- 
mogeneous mode (frnpoe ■ Therefore the normalization r 
does not introduce additional freedom into the solution; 
the homogeneous modes carry all the freedom. The value 
of ro will effectively be determined by determining the 
homogeneous modes, cither by specifying the material 
composition of the body to find its multipolc moments 
or by specifying global boundary conditions to find the 
pieces of the regular field. 

For the purposes of numerical calculations, one would 
likely use a normalization ro ~ 1Z (say, r$ = M in 
the case of an EMRI). Since the 'natural' normalization 
can be expected to be 2em, terms proportional to In e 
will likely appear in the homogeneous modes h^ Q Q ^ 

and h^iQ-^y corresponding to the logarithmic terms 

hirh^Q 1 \ and r hxr}v^\ 1 in the second-order field. 
This would mean lne terms would appear in the reg- 
ular field h^) '■ Whether or not that is the case can 
be determined analytically by evaluating the integrals in 
Eq. (|137p . But note that these potential lne terms, and 
any that appear at higher order in e, would not disrupt 
the ordering of the perturbations: e™(lne) 9 3> e™ (lne) 9 
for any n' > n, regardless of the values of q and q' . 

Note that Ref. 22] disallows In r terms in hf V s because 
— . W 

they will correspond to In e terms in the inner expansion, 

where lnr would become ln(ef). If the logarithms occur 
due to the standard retardation effects just described, 
then the lne terms appear instead in the outer expan- 
sion. Here, since I assume only the existence of a certain 
asymptotic series, rather than smoothness in e at e = 0, 
and since the logarithms are generic, I do not disallow 
them. A wave gauge in which they do not appear would 
certainly be advantageous, however, particularly at or- 
ders where logarithmic terms might lead to In e terms 



in a multipole moment I^ L , potentially obscuring the 
moment's meaning. 
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Appendix D: Second-order singular field 

The second-order singular field can be written as the sum of three pieces, = (_2) ^-(2) + (-i)^te) ^{2)^- Those 
three pieces are given below through order r, in Fermi- Walker coordinates with the mass dipole set to zero. I have also 
made the second-order field available in a Mathematica notebook at http:/ /www. personal. soton.ac.uk/ap8ell/second- 
order-ficlds.html. The notebook utilizes the tensor-manipulation package xAct. 

c _ 2) ftf 2) = (2e ab c S b a c , t - fS b B ab ) n a + §S°£ 6c n afec + r(jS a a b B ba - lS a a b e ba c a c ,t + 7S a a b e ac d a d , t n b c 

+ 2S a a b e bac a dit h cd + £S a a b B b c n ac + ^S a a b B a c n bc - §a b S b B ac n ac - %S b B bac n ac + \S b £ ac e bcd n a d 
+ IS b B acd h bacd + fS a a b B cd n bacd ) + 0{r 2 ), (Dl) 

(-2)^(2) = - 2e<ll3 r 2 Jni + ^^f^ + (\S c £ b d e a cd - f 5 S c £ ad e hcd ) n b + 3a a a b S d e bcd n c + \S b £ cd e bc ^ dl 

- S b £ cd e a b l h cdl + fa b a c S d e a d t h bcl + r(\S b e a bc a d , tt h cd + \S b e\ c a cM - ±S b B cd n a bcd - £S a B bc h bc 
+ §S b B b c h a c + §S b B ac h bc + ^S b B\ - \a b S c £ di eJh a bij + la b S c £ di e a Jn bdij + ^a a S b £ cd e bdl 

i 2 b Qccdi ~a , 3„bncc d ~ai ■ 13 _6 Qccdi a ~ , 69 _fc Qccad ~ i ■ 17 b qc cdi a » 

+ 2]-a 5 t e bc in d + ja b t b e cd in + j^a b t e C in bd + j^a b t e cd in b + bee bc n d i 

i 25_bccc da » i 7 b qc cad ~ i , 2 n b qc c da , 191 „6 Qccad r , 1 ob j ccdi^a 

+ g^a b tb e C in d — -^a b c e bc in d + j^a b c b e cd + -jg-a b c e bcd + ^b ej, c J £ n 

- i^e%^ c *n cd4J + i5 b e a b4 f c /n cd + %S b '£ a * e bdl n cd - fa b a c a d S l e\^n bcdj + fa a a b a c S d e cdl n b t 
+ 0(r 2 ), (D2) 

( _2)^ 2 fc ) = (2S d a^\ t e b \ d + ±S d 5 ab B cd - f S (a B b) c + fS c B ab ^j n c + ^S c B (a c n b ^ + 2 (V d S* (a + S c B d(a ^j h b \ d 

- 2S c B d h ab d - 2S c S ab B dl n cdl + r(fB {a c S b) a c + ^B c d a c S {a n b) d + fB cd a^ a S b} n cd ~ ^B d{a S b) a c h cd 



ff S^ a B b) cd h cd + ^B cdl S (a h b) cdl + ^B dl a c S {a n b) cd% - %B ab a c S c + §a c S c B d{a n b) d - ^B dl a c S c 5 



,,n 



f.b)d 



- fB ab c S c + l£ d{a e b) cd S c + ±B {a c a b) S c + ^£ dl S c ^ a c% h b \ - fB c d S c a (a h b ^ d + §S c B {a , 
+ IS c e cdt £ d ^ + §S c B d ^ a a b ^h cd + ^B ab d S c n c d - §£ d (-e b K t S c h/ - l£ di e cij 6 ab S c nJ + %B cdl 5 ab S c n dl 

- lB c dl S c h ab dl + p c B dl( - a n b \ dl + l£ di S c e^ cj n b ^ di - §B d « S c 6 ab n cdij ~ lB cd a c S ab S d ~ §B cd a c S d h ab 
+ £a c S d B( a d n b K + §S d a c B^ a c n b \ - § B ab a c S d n cd + fBj S d a c 6 ab h a - f 1 B c l a c 8 ab S d h dl 

- fB d l a c S d n ab cl + fa c S d B^ a n b \ dl - fB ij a c S ab , S d h aHj + 3a c S d a^\ t e b) dl h c l + ^ a e b \ d S c a d , t 

+ ^ a e b K d a ht S c n d ^+0(r 2 ), (D3) 

- &7TI / 

(-1)^(2) = \bm tt a a n a + r(^8m t \ tt + 2a b (5m* M + Sm tb a b ^ + ^5m tt a b a b - 8m ba £ ba + a a a b 5m ba 
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